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Abstract 


The  present  report  contains  further  analysis  of  the  radiation  data 
obtained  by  satellite  Explorer  IV,  1958  Epsilon,  especially  for  ARGUS 
charged  particle  geomagnetic  field  trapping.  The  present  work  was 
undertaken  to  supplement  the  previous  analysis  (Lundquist,  Naumann  and 
Weber,  J.  Geophys.  Res.  67.  4125  (1962)) for  the  naturally  trapped 
radiation. 

Directional  count-rates  (using  the  original  telemetry  data  from  the 
scintillation  counter,  Detector  A,  Channel  2)  as  a  function  of  time 
and  0  (angle  between  counter  and  plane  perpendicular  to  geomagnetic 
field)  are  plotted  and  analysed  with  comparison  of  Natural  and  ARGUS 
trapping  for  two  passes  of  the  satellite  over  Huntsville,  Alabama. 
Complete  analysis  of  all  data  acquired  by  Pass  4l4  (including  ARGUS 
Event  I  data)  is  presented  as  an  example  of  the  detailed  analytic  pro¬ 
cedure  which  is  required,  with  partial  analysis  of  Pass  4^4  (Event  II) 
included  for  comparison.  The  true  (corrected  for  deadtime,  effective 
counter  area  and  efficiency)  count-rate  vs  0  curves  were  reduced  ana¬ 
lytically  to  obtain  angular  distributions  of  the  directional  charged 

-2  -1  -1 

particle  flux  densities,  J(0‘)  (cm  sec  ster  )  vs  0'  (angle  between 
direction  of  J  and  counter  axis). 

The  results  indicate  that  the  Natural  angular  flux  density  distri¬ 
butions  are  roughly  constant  in  shape  for  nearly  the  same  region  in 
space  in  agreement  with  previously  published  results.  Por  ARGUS  re¬ 
gions  the  angular  distribution  of  the  directional  radiation  is  deter¬ 
mined  for  time  periods  in  which  the  background  omnidirectional  pene- 
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trating  radiation  (recorded  by  the  G-M  counters,  Channels  1  and  3)  is 
sensibly  constant.  Sudden  changes  in  the  omnidirectional  background 
counting  rate  are  shown  to  modify  the  shapes  of  the  count-rate  versus 
6  curves  which  correspond  to  individual  modulations  of  the  scintillation 
detector  count-rates  for  Pass  4i4  in  the  ARGUS  region.  It  is  concluded 
on  the  basis  of  the  limited  data  completely  analysed  that  the  corres¬ 
ponding  angular  distribution  of  the  directional  radiation  in  space,  and 
therefore  also  the  related  mirror-point  distributions,  for  ARGUS  trapped 
particles  are  affected  quite  markedly  by  the  omnidirectional  penetrating 
radiation  present  as  background.  ARGUS  particle  fluxes  observed  were 
1.0-1. 9  x  10®  and  .65  x  10®  cm"2sec"1ster“1  for  Passes  414  (Event  I)  and 
454  (Event  II)  respectively. 

Omnidirectional  count-rates  (Detectors  C,  unshielded  G-N  counter, 
Channel  3;  D,  shielded  G-M  counter,  Channel  1)  versus  time  are  presented 
for  Passes  414,  415,  4l6,  427,  440,  453,  454,  466,  and  479  following 
Event  I,  and  for  Passes  453,  454,  455,  466  and  479  following  Event  II. 

Some  broad  energy  dependence  of  the  omnidirectional  radiation  is 
shown  by  plots  of  both  the  ratios  of,  and  the  differences  between,  the 
unshielded  and  shielded  0-M  count-rates  which  were  obtained  for  the  same 
set  of  Explorer  IV  passes.  Time  decay  curves  of  the  omnidirectional 
ARGUS  trapped  radiation  were  constructed  for  constant  L  both  for  the 
unshielded  G-M  count -rates  and  the  ratio  of  unshielded  to  shielded 
0-M  count-rates.  The  ratio  curves  for  L  ■  I.72  and  1.73  (Pig.  3*4) 
indicate  an  unexplained  radiation  hardening-softening-hardening  sequence. 

The  Explorer  IV  roll-period  history  was  studied  in  some  detail  in  the 
region  of  the  reported  (Maumann,  1961)  anomalous  behavior.  Usiiqt  a 
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power  spectrum  analysis  of  the  count-rate  modulations  of  the  direc¬ 
tional  scintillation  counter  (Detector  A)  it  is  shown,  as  suggested 
by  J.  Bock  (1964),  that  there  is  no  anomalous  roll-period  disconti¬ 
nuity  if  beat  and  harmonic  frequencies  (of  the  roll  and  tumble  rates) 
are . recognized  and  eliminated  (see  Pig.  3.3  of  present  report). 

Considerably  more  data  remains  to  be  analysed  to  verify  more  posi¬ 
tively  or  modify  the  tentative  conclusions  concerning  directional 
particle  distributions  for  both  the  natural  (especially  pre  ARGUS) 
and  ARGUS  trapping,  to  determine  natural  and  ARGUS  mirror-point  dis¬ 
tributions  and  trapping  mechanisms,  and  to  explain  certain  features 
of  the  omnidirectional  radiation  energy  dependence. 


ill 


"Analysis  of  ARGUS/ Explorer  IV  Records" 

Chapter  1 

Introduction 

The  present  report  includes  a  further  analysis  of  the  large  amount 
of  geomagnetically  trapped  particle  radiation  (due  to  both  natural  and 
ARGUS  injection)  data  obtained  by  Explorer  IV,  satellite  1958  Epsilon, 
during  July  -  September  1958. 

A  previous  report  (Lundquist,  Naumann,  and  Weber,  J.  Geophys.  Res. 

67,  4125  (1962))  presented  an  analysis  of  the  naturally  trapped  radi¬ 
ation  as  detected  by  the  directional  scintillation  counter  (Detector  A, 
Channel  2)  on  board  the  satellite  and  telemetered  from  Explorer  IV  to 
the  Huntsville,  Alabama  ground  station. 

The  present  work  was  centered  on  ARGUS  type  trapping  and  was  in¬ 
tended  to  supplement  the  previous  report.  The  objectives  were  to  apply 
analysis  similar  to  that  employed  for  the  naturally  trapped  particles 
(op.  cltj  to  the  directional  ARCUS  trapped  particles  and  also  to  attempt 
further  analysis  of  the  omnidirectional  penetrating  radiations,  both  for 
natural  and  ARGUS  trapping,  using  the  data  obtained  by  the  G-M  detectors 
(Detectors  C  and  D,  Channels  3  and  1). 

The  following  chapters  of  the  present  report  develop  these  objectives 
and  present  all  results  so  far  obtained. 


Chapter  2 


Satellite  Explorer  IV  Instrumentation,  Motion  and  Telemetry 
2.1  Instrumentation.  Calibrations 


A  complete  description  of  the  detection  instrumentation  has  been 
published  by  Van  Allen  et  al,  (1959  a,  b).  A  block  diagram  showing 
the  arrangement  of  the  4  detectors  (A,  B,  C,  D)  and  the  corresponding 
telemetry  channels  2  and  5,  4,  3,  1  is  shown  by  Pig.  2.1.  A  sketch 
of  the  physical  arrangement  of  components  in  the  satellite  is  shown 
by  Pig.  2.2.  A  summary  of  the  characteristics  of  the  four  radiation 
detectors  is  given  in  the  literature  (Table  1,  Van  Allen  et  al,  1959b). 

Both  the  scintillation  counter  and  the  G-M  counter  count  rates  must 
be  corrected  for  deadtimes  (op.  cit.).  For  Detector  A,  the  plastic 
scintillator-photomultiplier  counter  and  associated  electronics,  the 
deadtime  correction  is  accomplished  by  the  equation 


c  *  Cobs  /  * 1  "  To  cobs  *  » 


(2.11) 


where  the  observed  counting  rate,  C  is  the  true 

counting  rate,  and  X  Is  the  deadtime  (92  microseconds). 
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For  Detector  C,  the  Anton  302  G-M  counter,  the  equation 


Cobs  *  c  6  T|  ° 


(2.12) 


where  r  (sore  properly  defined  as  deionization  time,  rather  than 
deadtime)  ■  62.5  ±  1.3  microseconds,  was  used  (op.  clt.).  A  plot  of 
Eq.  (2.12)  is  shown  by  Fig.  2.3.  A  discussion  of  the  actual  deadtimes 
used  in  the  present  analysis  follows. 
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In  the  present  work  both  the  T  and  T  values  of  Van  Allen 

0  I 

were  checked.  The  value  of  T  =  62.5  ±1.3  microseconds  was  con- 

*  4 

firmed.  However,  for  Detector  A  count-rates  near  a  maximum  of  10 

per  sec  were  observed  in  the  telemetry  records  (Pass  4i4).  Consider- 
4 

ing  10  per  sec  as  the  saturation  counting  rate  yields  (setting  de¬ 
nominator  of  Eq.  (2.11)  to  zero)  =  100  microsec.  In  turn, 

D 

Baicy  et  al  (1962)  determined  T  as  112  microsec.  Since  Baicy's 

0 

value  of  T  would  permit  a  maximum  Q  .  of  aboufc  9000  per 

0  ODS 

sec  only,  it  was  concluded  (for  the  present  work)  that  the  prototype 
detector  calibrated  by  Baicy  had  a  differing  from  that  of  the 

detector  actually  flown  in  Explorer  IV. 

It  was  decided,  in  the  present  analyses,  to  use  =  100  microsec 
for  Detector  A,  a  somewhat  arbitrary  decision.  Such  decision  is  how¬ 
ever  necessary  in  the  absence  of  detailed  calibration  information.  It 
should  be  emphasized  that  any  uncertainty  in  Tq  will  be  reflected 
in  the  shape  of  the  directional  flux  density  J(0‘)  curves  of 
Chapter  4.  The  maxima  of  these  J(®*>  curves,  near  « 0  • Min 
be  especially  sensitive  to  the  Tq  value  used. 

Additional  calibration  information  especially  for  the  directional 
scintillation  counter  (Detector  A)  was  determined  for  a  prototype  pack¬ 
age  (Baicy  et  al,  1962),  as  follows. 

1.  Detector  A.  The  pulse  height  discriminator  of  the  Detector  A  circuit 
(Pig.  2.1)  had  a  threshold  of  700-KeV  for  electrons.  The  deadtime 
for  the  circuit  was  measured  electronically  as  112  microseconds. 

The  discriminator  threshold  for  ganna-rays  was  between  0.66  and  1.17  MeV. 
The  ganwa  radiation  was  penetrating  enough  so  that  small  shielding  effects 
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were  not  observed.  No  multiple  scattering  into  the  detector  from 
other  parts  of  the  package  was  evident.  The  x-ray  calibration  es¬ 
tablished  a  1  MeV  photonic  threshold;  detector  window  response  curves 
demonstrated  the  existence  of  a  directional  plus  a  penetrating  omni¬ 
directional  component^  and  with  increasing  energies,  the  omnidirec¬ 
tional  component  obscured  the  directional  component.  The  calibration 
indicated  that  the  half-width  angular  response  curve  for  photonic  radi¬ 
ation  is  about  30#  broader  than  that  obtained  for  electrons  (Pig.  4.2). 

Over  the  measured  energy  range  (0.2  to  1.5  MeV),  the  angular  response 
curves  were  found  to  be  almost  independent  of  energy.  The  discrimi¬ 
nator  transmission  efficiency  €  and  effective  area,Agf  j  for  0-rays 
were 

€ » o,  €  A  ^||  *o  for  E<  700  KeV.-  (2.13) 

«•§.»%,  mm*  for  E>700  K6V  (2.14) 

The  calibration  data  noted  by  Eqs.  (2.13  and  2.14)  were  employed 
in  the  present  work  for  the  calculations  of  directional  electron  fluxes. 

For  an  average  photon  energy  of  1.25  MeV,  the  omnidirectional  factor 

(*  f)omnj  -  6.29  x  10-2  mm2  .  (215) 

2.  Detectors  C  and  D.  The  counters  were  identical  Anton  type  302  Q-M 
tubes  located  side  by  side  with  a  center  line  separation  of  3.6  cm, 
with  Detector  D  having  additional  shielding  (Van  Allen  et  al,  19590). 
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The  output  of  the  G-M  tubes  and  associated  circuitry  began  to  saturate 

at  approximately  2000  counts/sec  for  photonic  radiation.  An  analysis 

of  observed  counting  rate  for  pass  4l4  indicated  saturation  was  at 

approximately  58OO  counts/second.  A  further  increase  of 

radiation  flux  produced  a  decrease  of  count-rate.  The  G-M  tubes  were 

much  more  efficient  for  photonic  radiation  than  Detector  A.  No  counts 

were  obtained  for  electrons  up  to  1.5  MeV  due  to  the  satellite  shell 

stopping  power  of  1.2  gm/cm  .  Because  of  low  conversion  efficiency, 

it  was  not  possible  to  observe  any  bremsstrahlung  produced  by  incident 

electrons  on  the  satellite  shell.  The  omnidirectional  effective  areas 

for  Detectors  C  and  D  were  found  for  photons  (Baicy  et  al,  1962). 

2 

Detector  C  was  shielded  by  1.2  gm/cm  of  stainless  steel.  If  the 
shielding  were  aluminum,  the  extrapolated  range  would  be  given  as  a 
function  of  energy  by  Marshall  and  Ward  (1937) 

R«  0.526 E  “*  0.094,  (2.16) 

e 

where  is  measured  in  gm/cm  and  E  in  MeV.  Equation  (2.16)  yields 
the  absorber  thickness  for  10  percent  transmission  for  electrons  of 
energy  E.  The  range  for  50  percent  transmission  of  electrons  is,  same 
units  as  Eq.  (2.16), 

R  »  0.3«E  “  0*1  .  (2.17) 

1/2 

To  obtain  the  thickness  of  aluminum  equivalent  to  a  given  thickness  of 
another  absorber  of  nucleon  number  A,  the  rbsorber  thickness  is  multi¬ 
plied  by  the  factor  (Husain  and  Putnam,  1957) 


(2.18) 
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2 

The  1.2  gm/cm  of  stainless  steel  (  Z  =  26  )  therefore  are  equivalent 
2 

to  1.46  gm/cm  of  aluminum,  the  extrapolated  range  corresponding  to 
2.95  MeV  electrons.  For  'jO  percent  transmission  the  electron  energy 
is  4.1  MeV.  These  ranges  correspond  favorably  with  the  published 
values  (Table  1,  Van  Allen  et  al,  1959b). 

The  equilibrium  flux-energy  beta-ray  spectrum  of  U2^  was  reported 
by  Carter  et  al  (1959)  as 

N(E)  *  3.88  exp  (  -0.575E  -0.055E2  ),  (2. 

where  E  is  .the  0-ray  kinetic  energy  in  MeV,  and  N(E)  is  the  absolute 

differential  spectral  intensity  of  0-rays  in  electrons  per  fission  per 

MeV.  The  integrated  fluxes  have  been  obtained  by  numerical  integration 

(Van  Allen  et  al,  196?)  whence  the  fraction  of  electrons  of  energy 

_2 

greater  than  3  MeV  is  7.7  x  10  . 

2  2 
Detector  D  had  1.6  gn^cm  of  lead  in  addition  to  the  1.2  gVcm  of 

2 

stainless  steel.  The  total  shielding  is  equivalent  to  3.28  gm/cm  of 
aluminum  which  corresponds  to  an  extrapolated  range  of  6.2  MeV  electrons. 
The  fraction  of  electrons  of  energy  greater  than  6.2  MeV  is  1.6  x  10~^, 
approximately. 

The  efficiency  of  the  G-M  tubes  for  the  direct  counting  of  penetra¬ 
ting  electrons  and  for  counting  non-penetrating  electrons  via  inter¬ 
mediate  bremsstrahlung,  assuming  a  fission  spectrum  of  0-rays,  should  be 
Investigated.  Motz  and  Carter  (1963)  studied  the  response  of  an  Anton 
302  q-m  counter  to  a  fission  0-spectrum.  They  found  that  up  to  about 
2  govern  of  lead  absorber,  the  transmission  curves  were  as  expected  for 
0-transml88ion;  however,  for  greater  thicknesses  the  transmission  was 
enhanced  by  bremsstrahlung.  The  shielding  of  Detector  D  is  equivalent 
to  2.16  gs^cm2  of  lead  thus  indicating  bremsstrahlung  production  should 


be  considered. 
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2.2  Body  motion  of  Explorer  IV 


The  satellite  Explorer  IV  was  injected  into  orbit  July  26,  1958 
at  approximately  1500  UT  at  a  longitude  of  264.82  degrees  East,  a 
latitude  of'  12.95  degrees  North  and  an  altitude  of  368.2  km.  The 
initial  perigee  of  the  orbit  was  163  miles,  the  apogee  1380  miles, 
period  110.27  minutes;  the  inclination  to  the  equatorial  plane  50.29 
degrees.  Eccentricity  of  the  orbit  was  0.1297.  The  nodal  shift  per 
revolution  was  27. 9  degrees.  The  injection  pitch  flight  angle  was 
89*8  degrees,  the  yaw  flight  angle  0.8  degrees. 

Figure  2.2  is  a  sketch  of  the  satellite  showing  dimensions  and 
configuration.  The  total  weight  was  38.4  lbs  of  which  25.8  lbs  was 
instrumentation.  The  center  of  gravity  was  located  45  in.  from  the 
nose.  The  moment  of  inertia  about  the  long  or  Z1  -axis  was  0.50  in. 
lb  sec  r  about  the  X  -axis,  42.0  in.  lb  sec  .  The  cross  sectional 
area  was  30.7  sq  in.,  the  maximum  visible  projected  area,  458  sq  in. 

The  skin  of  the  satellite  was  stainless  steel  which  had  been  given  a 
fine  sandblast  resulting  in  an  emlsslvlty  0.54  and  0.46  for  the  wave¬ 
length  range  8-12  microns. 

Electrically  the  skin  of  the  satellite  was  divided  into  three  in¬ 
sulated  sections  as  shown  in  Fig.  2.2.  The  nose  cone  and  Instrument 
section  (designated  by  L  )  served  as  the  antenna  (asymmetrical  dipole) 
for  the  low  power  (10  mw)  108.00  me -sec transmitter.  The  two  middle 
sections  served  as  the  antenna  (symmetrical  dipole)  for  the  108.03 
mc-sec"1  high  power  (25  mw)  transmitter.  The  resulting  antenna  radi¬ 
ation  patterns  are  sketched  in  Pig.  2.4. 
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At  the  separation  of  the  second  stage  the  satellite  was  spinning 
about  the  z’ -axis  (Pig.  2.2)  at  a  rate  of  24  «  . radians/sec.  This 
Z‘-axis  spin  represents  a  maximum  energy  configuration.  As  rota¬ 
tional  energy  is  dissipated  (through  vibrational  friction  since 
the  satellite  is  not  completely  rigid)  conservation  of  angular  momen¬ 
tum  required  an  increase  in  moment  of  inertia  about  the  L-axis  to  com¬ 
pensate  for  decrease  in  angular  velocity.  Hence  the  motion  changes 
from  pure  spin  (or  roll)  about  the  (longitudinal)  symmetry  axis 
(A,  Pig.  2.2)  of  the  satellite  to  a  more  complicated  motion  of  spin 
and  tumble.  The  rate  at  which  this  transition  occurs  depends  on  the 
energy  dissipation  rate.  The  transition  from  the  initial  pure  spin 
to  nearly  pure  flat  or  propeller-like  motion  was  observed  to  occur  in 
about  9  days.  This  transition  was  confirmed  by  the  calculations  of 
Snoddy  (1939)  which  showed  that  the  temperature  variations  of  the 
satellite  could  be  explained  by  the  Increased  average  surface  area 
presented  to  the  sun  if  the  opening  angle,  between  the  pure  spin  axis 
and  the  spin-tumble  axis,  Pig.  2.5)  would  Increase  about  10  degrees 
per  day.  The  transition  from  maximum  to  minimum  rotational  energy  Is 
treated  by  Goldstein  (1939).  Slater  and  Prank  (1947).  and  others. 

Using  the  Euler  angles  as  described  by  Slater  and  Prank  (op.  clt.)  the 
following  important  results  can  be  derived 

L  *  lrwz,  /cos®; 

^*L/l  ,(ol i  )  ,  |  «  |  ; 

r  X  pox  r  9 

♦  ■  Ij./  Ix,  {Wf/COsB)  ,  (2.21) 


8 


In  Eqs.  (2.21)  the  primes  refer  to  the  satellite  fixed  rotating  coordi¬ 
nate  system  Pig.  2.5  and  the  angular  momentum  vector  is  considered 
to  be  constant  for  the  present  and  its  direction  used  to  define  a  space 
fixed  Z  -axis  as  shown  in  Pig.  2. 5.  where  x  and  Y  form  a  right  hand 
space  fixed  coordinate  system  with  Z,  If  energy  is  dissipated  the 

1  * 

Z  -axis  performs  a  torque-free  precession  about  L  with  precession  rate 
where  (f>  is  the  Euler  angle  measured  in  theX”Y  plane  (Pig.  2.5) 
and  is  the  angle  from  the  space  fixed  X-axis  to  the  line  of  nodes,  6  is 
the  angle  (  Z ,  Z*  )  or  the  opening  angle  of  the  precession.  The  third 
Euler  angle  \f/  is  the  angle  in  the  plane  of  rotation  about  Z*  from 
the  line  of  nodes  to  the  rotating  satellite  fixed  /  -axis.  By  the  first 
of  Eqs.  (2.21)  if  ^  .  decreases,  cos  Q  must  decrease  and  Q  increase. 

z  » 

The  precession  angular  velocity  <p  is  constant  and  may  be  evaluated 
from  the  second  of  Eqs.  (2.21),  the  initial  condition  when  rotation  is 
about  the  Z  -axis  only.  Since  the  ratio  of  the  transverse  moment  of 
inertia  to  the  longitudinal  is  84,  the  precession  frequency  is  1/84  of 
the  initial  spin  rate  of  24  it  radlans/sec.  or  0.286  it  radlans/sec. 

(period,  7  secs).  When  the  motion  has  degraded  to  the  pure  propeller- 
type  (tensed  "tumble"  as  opposed  to  "spin-tumble")  motion,  the  corres- 

m 

ponding  frequency  is  tensed  the  tumble  frequency  (  <f>  is  the  corres¬ 
ponding  tumble  angular  velocity)  and  has  been  measured  as  described  in 
Sec.  3.3  by  observing  the  frequency  of  r-f  signal  fading  at  such  times 
when  the  satellite  la  properly  orientated  relative  to  the  receiving 

station.  The  satellite  continues  to  have  a  slight  angular  spin  velocity 
* 

about  the  Z*  -axis.  The  period  of  this  spin  or  roll  motion  was  ob¬ 
served  to  vary  in  magnitude  during  the  lifetime  of  the  satellite  from  1 
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to  approximately  100  secs.  The  combination  of  the  spin  and  tumble 
motions  result  in  the  sweeping  of  the  onboard  radiation  counters  over 
almost  all  directions  in  space;  in  particular  the  directional  counter 
(Detector  A)  was  able  to  observe  counting  rates  at  all  angles  relative 
to  the  geomagnetic  field  direction  between  0  and  rc/2.  The  analysis  of 
directional  radiation  flux  in  space  is  seen  therefore  to  depend  on  ob¬ 
taining  with  good  precision  the  orientation  of  the  scintillation  counter 
with  respect  to  the  magnetic  flux  density  ft  as  a  function  of  time. 
During  the. time  of  acquisition  of  the  data  analyzed  in  the  present  report 
the  satellite  remained  in  the  flat  propeller-like  configuration  and 
rotation  about  the  angular  momentum  vector  will  always  be  referred  to 
as  tumble  (velocity  or  rate,  )  and  motion  about  the  Z*  -axis  will 
be  referred  to  as  spin  (or  roll)  (velocity^or  rate,  ^  ),  as  has  been 
noted  above. 
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2.3  Interpretation  of  satellite  Explorer  IV  telemetry  data 

Figure  2.6  is  a  reproduction  of  a  portion  of  a  record  of  the  tele¬ 
metry  information  received  from  the  satellite  at  Huntsville,  Alabama 
on  August  27,  1958  following  04i6  UT,  for  the  region  of  ARGUS  trapping. 
The  radiation  counter  data  shown  was  recorded  on  magnetic  tape  and 
played  back  for  the  pen  on  paper  recording  using  electronic  filters 
suitable  for  reducing  noise.  The  variations  in  the  signal  represent 
the  abrupt  changes  in  frequency  of  the  audio  channel  modulating  the 
carrier  of  one  of  the  transmitters.  For  example,  the  first  line  of 
information  (Channel  1)  represents  audio  switching  back  and  forth 
from  370  to  430  cycles/sec.  Each  frequency  switch  follows  the  accu¬ 
mulation  of  32  counts  from  the  shielded  G-M  counter,  Detector  D 
(reference.  Fig.  2.1).  Similarly,  changes  in  the  second  line  (Channel  2) 
of  information  follow  the  accumulation  of  1024  counts  (per  frequency 
switch)  of  the  scintillation  counter,  Detector  A;  the  third  line 
(Channel  3)  shows  frequency  switching  following  the  accumulation  of 
1024  counts  (per  switch)  for  the  unshielded  G-M  counter,  Detector  C. 
Detector  A  Is  the  directional  counter  and  the  rapid  changes,  or  modu¬ 
lations,  of  counting  rate  of  Channel  2  correspond  to  the  changing 
orientation  of  the  counter  axis  with  respect  to  the  magnetic  field 
direction  as  noted  In  Sec.  3.1. 

Figure  2.7  illustrates  the  r-f  field  strength  lnfonnation  received 
at  the  Huntsville  ground  station.  This  information  is  used  to  determine 
the  orientation  of  the  satellite  and  Its  radiation  counters  relative  to 
the  geomagnetic  field  direction  (see  Chap.  3).  Each  telemetry  record 
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also  contains  a  recording  (not  shown  in  Pigs.  2.6,  2.7)  of  the  time 
signal  from  WWV,  permitting  real  time  determinations  to  0.05  sec. 

Records  similar  to  Figs.  2.6  and  2.7  were  reduced  for  Passes  4i4, 


415,  4l6,  427,  428,  440,  44l,  453,  454,  455,  466,  479,  492,  505,  and 
518  by  using  a  Telecomputing  digitizer.  An  initial  reference  time  t 

c 

was  established  for  each  time  interval  T  (24  secs).  N  is  the  number 
of  digitizer  counts  per  24  sec  interval  and  (  H  |  “  n  j  )  is  the 
number  of  digitizer  counts  per  frequency  change.  These  data  were 
supplied  to  an  IBM  1620  computer  along  with  the  scaling  factor  S 
(noted  on  Pig.  2.1)  and  used  to  compute  the  observed  count  rate 
(counts/sec),  as  follows 


C0bt  *  NS/2T(a^|”  n. ), 

and  the  corresponding  time  t^  is  calculated  by 

ti',o+(ni  +  "i+l*  (T)/2N 

S  is  divided  by  2  since  S  is  the  number  of  counts  accumulated  during 
a  full  cycle  of  two  frequency  switches.  The  Fortran  program  used  for 
the  calculations  of  Eqs.  (2.31)  and  (2.32)  is  designated  as  Program  2.1 
(all  computer  programs  used  in  the  present  work  are  available  at 
Physics  department.  Saint  Louis  University,  St.  Louis,  Missouri).  The 
count  rate  vs  time  computer  output  for  Program  2.1  was  plotted  with  a 
Dymec  data  plotting  system.  Figure  2.9  Includes  the  resulting  observed 
count-rate  vs  time  for  the  directional  scintillation  counter  (Channel  2, 


(2.31) 


(2.32) 
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Detector  A)  as  Explorer  IV  traversed  a  region  of  naturally  Injected 
geomagnetic  ally  trapped  particles  a  short  time  (about  6  minutes)  be¬ 
fore  the  satellite  entered  the  artifically  injected  ARGUS  shell 
(Pass  4l4).  Figures  2.10  and  2.11  are  similar  plots  for  this  ARGUS 
penetration  (telemetry  playback  record  shown  by  Fig.  2.6). 

Figures  2.9  and  2.10  include  also  (Upper  curves)  plots  of  COS^$ 

(  9  >  angle  between  scintillation  counter  axis  and  the  plane  per¬ 

pendicular  to  the  geomagnetic  field  direction)  vs  time  following  the 
suggestion,  of  J.  Bock  of  the  Ballistics  Research  Laboratories;  9 
as  a  function  of  time  is  determined  by  Eq.  (3.11).  It  is  interesting 
to  note  that  the  modulations  of  the  observed  count  rates  vs  time 
(Figs.  2.9,  2.  10  Lower)  are  rather  well  represented  by  the  trigono¬ 
metric  function  (Figs.  2.9.  2.10  Upper)  since  the  calculation  of 
by  Eq.  (3.11)  is  based  on  experimentally  determined  parameters. 

As  noted  In  the  caption  for  Fig.  2.10,  the  Detector  A  count-rate 
modulations  were  individually  (a  through  q,  Figs.  2.10  (Lower)  and 
2.11)  smoothed  and  least-square  fitted  by  4th  degree  polynomials, 
using  computer  Program  2.2,  of  the  form 

A  ♦  A  t  ♦  ...  At"  <*•») 

o  i  n 

To  Interpret  the  nulls  in  the  ground  station  r-f  field  strength 
traces  shown  by  Fig.  2.7  it  is  helpful  to  consider  two  cases  with 
reference  to  the  satellite  antenna  radiation  patterns  shown  by  Fig.  2.4. 
First,  consider  the  receiver  antenna  in  the  plane  of  tumble  of  the 

A  A 

satellite  with  Uj.  L  ®s  shown  by  Fig.  2.8  Upper.  Hulls  in  the. 
field  strength  pattern  for  each  transmitter  occur  whenever  the  nose  or 
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tail  of  the  satellite  points  toward  the  receiver  producing  the  field 
strength  pattern  shown  by  traces  6  (due  to  108.03  mc-sec’1  transmitter; 
reference,  Fig.  2.4  Lower)  and  7  (due  to  108.00  mc-sec  1  transmitter; 
reference,  Fig.  2.4  Upper) .  Second,  consider  the  receiver  antenna  in 
a  plane  perpendicular  to  the  plane  of  tumble  of  the  satellite  with 

A  A 

U  ||  L  (Fig*  2.8  Lower) .  The  satellite  antenna  radiation  patterns 
shown  by  Fig.  2.4  produce  no  nulls  for  either  transmitter  antenna.  In¬ 
stead,  because  both  the  transmitter  antennas  and  the  receiving  antenna 
are  linearly  polarized,  nulls  occur  twice  for  every  revolution  of  the 

A 

satellite  about  L  when  the  plane  of  polarization  of  the  transmitting 
antennas  and  the  receiving  antenna  are  perpendicular.  These  polarization 
nulls  are  characterized  by  lack  of  the  extra  side  nulls  shown  for  the 
"nose"  nulls  by  Fig.  2.7  Lower.  Polarization  nulls  are  not  seen  in 

A 

Pig.  2.7  because  the  line-of-sight  vector  U  was  perpendicular  approxi- 

A 

mately  to  the  angular  momentum  vector  L. 

Naumann  (1961)  was  able  to  use  r-f  field  strength  information  as 
described  in  the  preceeding  paragraph  together  with  known  directions 

A 

for  the  line-of-sight  vector  U  (refer  Eq.  3.28)  from  the  receiving 
station  to  determine  the  time  history  orientation  of  the  angular  mo- 

A 

mentum  vector  L.  On  August  18,  1938  the  antenna  system  at  the 
Huntsville  ground  station  for  the  108.03  me  receiver  was  modified  to 
font  a  variably  polarized  array  which  changed  its  polarization  as  the 
satellite  passed  overhead.  This  arrangement  allowed  use  of  the  vari¬ 
ation  of  the  amplitude  of  the  nulls  with  the  variation  of  the  plane 
of  polarization  of  the  receiving  antenna  to  yield  additional  and  more 

A  A 

accurate  information  on  the  direction  of  L.  The  directions  of  L  used 
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to  analyze  the  passes  in  this  report  are  those  of  Naumann  (1961). 

The  direction  of  L  was  found  to  be  not  constant  as  expected  but  varied 
slowly  from  day  to  day;  however  for  telemetry  acquisition  periods  up 
to  some  20  min  the  direction  could  be  considered  constant.  Naumann 
(1962)  has  suggested  an  explanation  of  the  observed  variation  of 
angular  momentum  in  terms  of  the  interaction  of  the  magnetic  moment 
of  the  satellite  with  the  earth's  magnetic  field.  The  orientations 

A 

of  L  determined  in  this  way  checked  with  Snoddy's  (1959)  prediction 
of  variations  in  temperature  due  to  the  variation  in  satellite  pro¬ 
jected  area  toward  the  sun  and  orbital  trajectory  variations  due  to 
presentation  of  different  frontal  drag  orientations. 

The  rate  at  which  r-f  nulls  occurred  when  the  line-of-sight  vec- 

*  A 

tor  U  was  nearly  in  the  plane  of  tumble  variation  permits  the  cal¬ 
culation  of  the  satellite  tumble  rate  if  the  variation  due  to  motion 
along  the  trajectory  is  corrected  for  (see  Sec.  3.3.  Eq.  3. 32).  The 
tumble  period  actually  observed  is  6. 90  sec  approximately  (see  Pig. 
2.7),  nearly  7  sec  as  calculated  (Eqs.  2.21,  Sec.  2.2). 


15 


Pig.  2.1  Block  diagram  of  radiation  Detectors  A,  B,  C,  D  and 
corresponding  telemetry  Channels  2  &  5»  4,  3,  1.  The 
numbers  2048,  16,  DC,  2048,  64  designate  the  scaling  factors 
for  a  complete  cycle  (2  switches)  for  each  channel. 

(Van  Allen  et  al,  1959) - 


lo 


HIGH  ROVER 
25  NV  AM 

TRAMSMITTER 

108  03  MC 


Fig.  2.2  Sketch  showing  shape,  size  and  packaging  of  the  Explorer  IV 
satellite. 

Upper.  Complete  satellite  in  3  sections:  nose  cone  and 
instrumentation  section  (together  designated  by  L) ,  propellant 
and  exhaust  nozzle  (remaining  length). 

A 

Lower.  Nose  cone  and  instrumentation  section  only.  A, 
satellite  symmetry-axis  vector  (in  direction  of  Z  );  P, 

Detector  A  (directional  scintillation  counter)  axis  direction. 


17 


•0.39  in. 


ELECTRONICS,  BATTERIES 
&  HIGH  POWER  TRANSMITTER 


Fig.  2.3 

c 


Plot  of  Eq.  (2.12),  true  omnidirectional  counting  rate 


vs  observed  counting  rate  Cq|j5 


for  Detector  C, 


Anton  302  G-M  counter. 
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RUE  COUNTING  RATE  C  (SEC 


Pig.  2.4.  Antenna  radiation  patterns  for  transmitters  of  Explorer  IV 

(Naumann,  1961). 

Upper.  Pattern  of  asymmetrical  dipole  for  108.00  mc-sec”* 
low  power  (10  mw)  transmitter. 

Lower.  Pattern  of  symmetrical  dipole  for  108.03  mc-sec’^ 
high  power  (25  mw)  transmitter. 

Both  patterns  are  symmetrical  about  the  long  symmetry  axis 
of  the  satellite. 
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V 


Pig.  2.5  Important  geometry  for  describing  satellite  spin  (or  roll) 
and  tumble.  VJ  are  fixed  relative  to  satellite  (see 

Pig.  2.2),  with  t  as  spin  axis.  XYZ  are  space  fixed  with  Z 
(spin-tumble  axis)  taken  as  the  axis  of  (constant)  angular  momen¬ 
tum  L.  are  Euler  angles. 
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Pig.  2.6  Portion  of  ,:playback"  record  of  the  telemetry  tape, 
received  (August  27,  1958;  04i6  UT)  at  Huntsville,  Alabama 
for  Pass  4l4  of  Explorer  IV  Satellite. 

The  decreased  amplitude  of  the  recording  for  the 
middle  portion  of  Channel  1  Indicated  saturation  of  the 
shielded  G-M  counter  by  the  trapped  ARGUS  radiation. 

The  total  number  of  telemetry  channels  recorded  for 
records  such  as  Pig.  2.6  are  the  following  (3  channels  only 
appear  in  Pig.  2.6) 

Channel  1  - 
Channel  2  - 
Channel  3  - 
Channel  4  - 
Channel  5  - 
Channel  6  - 
Channel  7  - 
Channel  8  - 


Detector  D  (reference.  Pig.  2.1) 

Detector  A  (reference.  Pig.  2.1) 

Detector  C  (reference.  Pig.  2.1) 

Detector  B  (reference,  Pig.  2.1) 

Detector  A  (reference.  Pig.  2.1) 
r-f  pattern,  108.03  me -sec"1  transmitter 
r-f  pattern,  108.00  me -sec-1 
time  (WWV  time  signals). 
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Pig.  2.7  Sketch  (smoothed) of  r-f  field  strength  modulations  of 
the  Huntsville  ground  station  receiver  AGC  showing  variations 
due  to  rotation  of  the  satellite  antennas  (and  of  the  satellite 
Itself)  relative  to  the  line-of-sight  from  receiver  to  satellite. 
The  patterns  correspond  to  the  alignment  shown  by  Pig.  2.8  Upper; 
namely. 

Trace  6  is  for  the  108.03  me -sec "high  power"  transmitter, 
radiation  pattern  shown  by  Pig.  2.4  Lower ;  Trace  7  for  the  108.00 
mc-sec”1  "low  power"  transmitter,  radiation  pattern  Pig.  2.4 
Upper. 

N,  T  designate  "nose"  and  "tail"  ends  of  the  satellite 
respectively. 
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TIME 


Pig.  2.8  Sketch  illustrating  alignment  of  Explorer  I V  Satellite 
relative  to  ground  station  (Huntsville,  Alabama)  receiver 
antenna,  RA. 

A  A 

Upper.  U  L  L  Linearly  polarized  receiver  antenna 

A 

in  the  plane  of  tumble  motion  of  satellite.  U  ,  unit  vector 
defining  receiver  station  -  to  -  satellite  line-of-sight,  in 
plane  of  tumble. 

A  A 

Lower.  U  ||  L  Linearly  polarized  receiver  antenna 

A 

with  U  directed  perpendicular  to  plane  of  tumble. 

(A,  satellite  symmetry  axis  vector  direction;  reference 
Pig.  2.2). 
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Pig.  2.9  Lower,  observed  directional  count-rate  vs  time  (scintil¬ 


lation  counter  Detector  A,  Channel  2,  Pass  4l4)  for  naturally 

injected  trapped  particle  radiation  compared  with 

Upper,  variation  of  the  arbitrarily  chosen  trigonometric 
4 

function  cos  9  (6,  angle  between  scintillation  counter  axis 
and  the  plane  perpendicular  to  the  geomagnetic  field  direction) 
vs  time;  9  is  calculated  by  Eq.  (3.11)  and  such  calculations 
are  based  on  experimentally  determined  parameters. 
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UT,  August  2  7  (  mi  n  ) 


Pass  414)  Tor  first  half-minute  of  ARGUS  shell  penetration.  Crosses,  observed  counting  rates;  solid 


V°° 


( t_oas)sqoo 


August  27  (min) 


Fig.  2.11  Observed  directional  count-rate  vs  time  (scintillation 
counter,  Detector  A,  Channel  2,  Pass  414)  for  the  second 
half -minute  of  ARGUS  shell  penetration. 

Figure  2.11  follows  Fig.  2.10  immediately  in  time. 

The  descriptive  details  of  Fig.  2.10  caption  apply  also 
to  Fig.  2.11.  The  variations  in  the  omnidirectional  "pene¬ 
trating"  background  radiation,  solid  curve  B,  are  to  be  noted 
especially. 
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UT  August  27  (min) 


Chapter  3 


Determination  of  Count-Rate  Co^s  as  Function  of  the  Angle  6  Made  by 
the  Scintillation  Counter  with  the  Plane  Perpendicular  to  the  Geo¬ 
magnetic  Field 

3.1  Satellite  coordinate  systems  and  geometry 

The  directional  particle  flux  at  any  satellite  position  on  a  geo¬ 
magnetic  field  line  is  related  simply  to  the  mirror  point  distribu¬ 
tion  (Weber  et  al,  1962).  The  first  step  in  the  calculation  of  the 
directional  flux  at  a  point  in  space  is  the  measurement  of  the  count- 

A 

ing  rate  Cobs(6)  as  a  function  of  the  angle  6  between  P  ,  a  unit 
vector  along  the  axis  of  the  directional  counter,  and  a  plane  perpen¬ 
dicular  to  the  geomagnetic  flux  density  B  (unit  vector).  By  Fig.  4.1, 

6 -orcsintPB),  or  6  *arcos  (P'tfj, 

*  A 

where  BJ.  Is  a  unit  vector  perpendicular  to  B;  and,  for  pure  tumble 
motion  of  the  satellite 

6  *  orctin  Cb|  sin(^t+8)-  bf  cos(^t  +  8)  sin(^t  +  y)+ 1>9  co*^t+8)cos(^t+  y)3 

(3.11) 

In  Eq.  (3.11),  ir/Z-6t  the  tumble  rate  $  and  the  spin  or 

roll  rate  ^  are  defined  In  Sec.  2.2,  Y  and  6  are  Initial  phase 
angles  defined  by  Fig.  3.1*  and  b^,  b2,  bj  are  the  components  of 
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vector  B  on  the  X  ,  Y  ,  Z  axes  of  the  satellite  inertial  coordi- 
s  s  s 

nate  system,  Fig.  3.1.  From  Eq.  (3.11),  6'  can  be  calculated  as  a 

function  of  time  and  compared  with  the  counting  rate  cobs(6)  as  has 

been  discussed  in  Sec.  2.3,  and  shown  by  Figs.  2.9  and  2.10,  pro- 
•  >  1  • 

vided  <J>  ,  are  known  and  the  phase  angles  Y  and  6  evaluated 

at  some  reference  time  tQ. 

The  geometry  of  the  present  analysis  involves  several  coordinate 
systems.  Equation  (3.11)  is  calculated  in  a  Satellite  Inertial  Coor¬ 
dinate  System  (Fig.  3.1),  in  which  X  has  the  direction  of  the 
angular  momentum  vector  L  (pure  tumble  axis  also)  considered  constant 
during  a  some  20-minute  interval  of  telemetry  data  acquisition;  the 
Z  direction  is  defined  by 

2t-  L*B./S,  where  s*  o.izi 

A 

and  Bq  Is  a  unit  vector  in  the  direction  of  the  geomagnetic  flux  den¬ 
sity  at  time  tA;  and  v_  is  determined  by 

U  S 


\  *  [tL*  B0«  L]  /S  .(E^-LK)/3,  where 


A  A 

K*  L  B0 


(3.13) 


The  Satellite  Inertial  Coordinate  System  of  Pig.  3.1  is  determined 

A 

such  that  B  is  in  the'  X  -Y  plane  at  t*t  . 

s  8  o 

The  axes  XA ,  Y  Z  of  the  Satellite  Inertial  Coordinate  System 

S  8 1  8 

are  to  be  unit  vectors  in  the  Vernal  Equinox  Coordinate  System  (Sec.  3.2) 
which  is  the  basic  reference  coordinate  system  used  in  the  present  work. 
In  the  Vernal  Equinox  Coordinate  System  the  X  axis  points  to  the  Vernal 
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Equinox  of  1951 »  the  y-  axis  is  in  the  earth's  equatorial  plane  and 
the  axis  coincides  with  the  axis  of  rotation  of  the  earth  at 
equinox.  The  components  b^,  b^  of  Eq.  (5.11)  in  the  Satellite 
Inertial  Coordinate  System  may  be  calculated  at  any  time  by 

b,*  B*„ 

be*  B-Yt,  (3.i4) 


it  being  noted  that  b^  is  zero  at  tQ. 

Figure  3.1  indicates  also  that  the  tumble  phase  angle  Y  is  defined 

by  the  direction  of  the  satellite  symmetry  vector  A,  (see  also  Fig.  2.2). 

A  third  coordinate  system,  a  Satellite  Fixed  Coordinate  System 

(Fig.  3.1)  is  used  to  define  6  ,  the  spin  or  roll  phase  angle.  In 

this  system  is  space-fixed  in  the  direction  of  L,  Yp  is  in 

the  direction  of  vector  A  (which  rotates  in  a  plane  perpendicular  to 

L  for  the  pure  tumble,  propeller- like  motion)#  Angle  8  is  the  angle 

between  P  and  Z  at  tine  t_,  L  is  supplied  by  Naumann  (1961)  in 
r  u 

A 

the  Vernal  Equinox  Coordinate  System,  B  can  be  converted  from  the 
ephemerls  system  of  coordinates  to  the  VEC  System  (Sec.  3.2)  and 
the  Satellite  Inertial  Coordinate  System  is  established  at  an 
initial  reference  t^  time  for  each  pass.  The  direction  of  the 

A  A 

vector  A  is  determined  from  the  direction  of  U  the  llne-of-sight 
vector  (see  Eq.  3.28)  from  the  ground  station  to  the  satellite  at 
the  tine  tQ  of  an  r-f  null  by 

A*  ♦  Bl*U)*L]/|LxU| 

where  the  +  sign  is  used  for  a  tall  null,  minus  sign  for  nose  null. 
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The  phase  angle  Y  is'  then  the  arccosine  of  the  dot  product  of  A  and 

Y  .  V/ith  the  various  parameters  of  Eq.  (3.11)  determined  e  may  be 

s 

calculated  for  any  time  following  t^  and  compared  with  the  count  rate 
C  .  evaluated  from  the  polynomial  expressions  as  discussed  in 

0  DS 

Sec.  2.3. 
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3.2  Conversion  of  1958  Epsilon  orbital  data  to  the  Vernal  Equinox 
Coordinate  System 

The  final  reference  coordinate  system  used  in  the  present  work  is 
the  Vernal  Equinox  Coordinate  System  (VE)  defined  and  briefly  discussed 
in  the  preceeding  Sec.  3.1*  Hence  it  is  necessary  to  convert  the  geo¬ 
magnetic  field  information  of  the  1958  Epsilon  Orbital  Data  Series 
(Smithsonian  Astrophysical  Observatory,  1959)  to  the  VE  Coordinate 
System.  Table  3.1  is  a  sample  listing  of  1958  Epsilon  Orbital  Data 
Series  (op.  cit.,  SAO,  1959)  for  Explorer  IV  Pass  4l4. 

In  order  to  compute  the  geomagnetic  field  components  &y  and  B 

in  VE  coordinates.  Program  3.1  was  used  to  convert  the  orbital  infor¬ 
mation  of  Table  3.1  to  the  form  shown  in  Table  3.2  where  time  is  ex¬ 
pressed  in  minutes  UT,  latitude  is  geocentric  latitude,  and  B  ,  B  r 

x  y 

and  B_  are  the  components  of  the  B  unit  vector  in  the  VE  Coordinate 

a 

System.  Program  3*1  yields  the  converted  data  at  1  min  time  Intervals 
and  also  computes  6th-order  polynomial  fits  to  these  data  (latitude, 
longitude,  altitude,  VE  geomagnetic  field  components  and  satellite 
velocity  components)  so  as  to  provide  interpolation  for  any  time  within 
the  30  min  interval  fitted  by  the  polynomials. 

Some  details  of  the  conversion  of  the  Explorer  IV  orbital  data  to 
Vernal  Equinox  coordinates  follow.  Calculation  of  the  vector  ?. 

(Pig.  3.2)  from  the  center  of  the  earth  to  the  Huntsville  ground  station 
and  the  vector  r  from  the  center  of  the  earth  to  the  satellite  In 
Vernal  Equinox  coordinates  Is  included  also. 

The  satellite  orbital  data  (see  Table  3.1)  lists  the  geographic  lati¬ 
tude  and  longitude  of  the  intersection  point  of  the  international 
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ellipsoid  and  the  radius' vector  r  to  the  satellite.  The  height 
listed  in  Table  3.1  is  the  distance  from  this  intersection  point  to 
the  satellite.  The  geomagnetic  flux  density  Table  3.1  is  in  rectan¬ 
gular  topocentric  coordinates  relative  to  a  spherical  earth's  surface 
with  north  N,  east  E  and  down  D  components  calculated  from  the  expan¬ 
sion  of  the  earth's  magnetic  potential  in  spherical  harmonics  using 
Finch  and  Leaton  coefficients  of  1955* 

The  latitude  most  commonly  listed  on  maps  and  in  tables  is  called 
the  geodetic  latitude  A  (Fig.  3.2).  This  is  the  angle  between  the 
perpendicular  to  the  tangent  to  the  international  ellipsoid  (Fig.  3*2) 
and  the  plane  of  the  equator.  Table  3.1  lists  the  geographic  latitude, 
defined  by  the  angle  between  a  local  plumb  bob  and  the  plane  of  the 
equator,  includes  the  effects  of  local  gravitational  anomalies  and  the 
local  centrifugal  force.  The  difference  between  geographic  and  geo¬ 
detic  latitude  is  small  and  called  the  station  error.  For  all  calcu¬ 
lations  in  this  report,  geodetic  and  geographic  latitude  A  are  taken 
to  be  the  same.  For  purposes  of  vector  calculation  and  conversion  to 
the  VE  Coordinate  System  it  is  necessary  to  use  the  geocentric  latitude 
A*  (Pig.  3.2),  which  is  related  to  geodetic  latitude  by 


ton  A1*  tonA  b2/a*;  b2/a2»  l/l.0043, 


(3.21) 


and  a  and  b  are  the  equatorial  and  polar  radii  respectively  of  the 
international  ellipsoid.  The  distance  r  of  the  satellite  from  the 
center  of  the  earth  in  earth  radii  is  equal  to 

r*  b/(a*smf  A'i-b*  cos^y4,  height  I958«/  o  , 


(3.22) 
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where  height  ^53  £  is.  altitude  listed  in  1958  Epsilon  Orbital 

Data  Series  (SAO,  1959)*  R  the  radius  vector  from  the  center  of 
the  earth  to  a  station  in  VE  coordinates  in  terms  of  the  geodetic 
latitude  A  and  the  longitude  ©  is  given  by 

Rx=(C+H)cosA  cos© 

Ry  =  (C+H)cosA  sin©  (3.23) 

Rz  *(S+H)  sin  A, 

where  longitude  angle  ©  is  the  hour  angle  to  the  Vernal  Equinox 
given  by 

0  *  0  j  +  0.0043752695  (t-tj)  +  ©E , 

in  which  @  •  is  the  angle  between  the  Greenwich  meridian  at  t.  and 

Ql  1 

t  is  in  minutes  after  ti#  ©^  is  the  east  longitude  of  the  station 
point;  and 

C*  ^l-(2f-fe)sinA)^*,  0.25) 

H*  receiver  antenna  height  above  the  international  ellipsoid  in 

earth  radii  units  f  ■  <n-b)/a  *  I  /297.0,  and  cO-f)*S. 


The  vector  r  from  the  center  of  the  earth  to  the  satellite  is  given 
in  VE  coordinates  by 

rx  •  h  cos  cos  6t 

r  •  h  cos  A,  sin  0S  (3.26) 

rt  •  h  sinA', 
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is  the 


B 


VE 


where  Ag  is  the  geocentric  latitude  of  the  satellite, 
east  longitude  of  the  satellite  from  the  Vernal  Equinox.  Thus,  the 

-V 

satellite  vector  r  is  calculated  using  A'  from  the  converted 
ephemeris  while  the  station  vector  R  is  calculated  using  A  in 
Eq.  (3.23)  which  Includes  correction  for  the  station  height  above  the 
geoid  and  the  difference  between  geodetic  and  geocentric  latitude. 

The  geomagnetic  field  strength  is  converted  by  Program  3.1  to  the 
VE  system  using  the  matrix  equation 

/BA  /I  0  0\  /cosA'  0  -sm^x/bA 

=  |byW°  cos8t  -sin0,  J  j  0  I  o  Jl  bf  V3'271 

\BJ  \°  cosSj  \  sin  A‘  0  cos  A'  J\bJ 


where  bjj,  bp  and  bp  are  the  normalized  (divided  by  B  magnitude)  of 
i  Bg  •  and  B  Once  B  is  determined  in  Vernal  Equinox  coordi¬ 
nate,  Eqs.  (3.12),  (3.13),  (3.14)  may  be  applied  to  calculate  b^,  b2 

and  b,  for  use  in  Eq.  (3.11). 

3 

A 

The  ground  station  to  satellite  line-of-sight  vector  U  from  which 

A 

the  satellite  symmetry  vector  A  is  calculated,  Eq.  (3. 15) •  is  given  by 


(3.28) 


It  remains  to  evaluate  the  tumble  and  spin  or  roll  rates  t  and  <1* 
and  the  tumble  and  spin  or  roll  phase  angles  y  and  6  .  The  deter¬ 
mination  of  these  parameters  is  considered  in  Sec.  3.3. 
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3.3  Determination  of  tumble  and  spin  or  roll  rates  and  phase  angles; 

•  • 

Y,  »•  The  spin  or  roll  period  anomaly.  Closure  — 

search. 

The  tumble  angular  velocity,  or  rate,  may  be  determined  as 

suggested  by  Naumann  (1961)  from  the  observations  of  the  period  of 
r-f  nulls  (Sec.  2.3)  correcting  for  the  orbital  motion  of  the 
satellite  as  is  indicated  by  Eq.  (3.32)  below.  The  direction  of  A, 
Pig.  2.2  may  be  determined  by 


where  and  0^  are  A  and  U  respectively  at  the  time  of  the  1th  null. 
The  angle  4  through  which  the  satellite  has  tumbled  between  the 
first  and  the  ith  null  is 

«  A  A 

*  2nir  -  arcos  (A,*  Aj) 

where  n  is  the  number  of  tumble  revolutions  between  the  first  and  the 
1th  null.  The  sign  of  the  arccosine  tera  depends  on  the  sense  of  the 
tumble  rate  and  the  direction  of  orbital  angular  momentum  vector. 
Equation  (3.32)  may  be  used  to  evaluate  b  .  However,  for  the  present 
report  the  tumble  rate  as  determined  by  Bock  (19^4)  was  used.  Bock  was 
able  to  determine  the  tumble  rate  to  5  decimal  place  accuracy  by  using 
an  analog  method.  For  Pass  414  Bock  (op.  cit.)  evaluated  the  tumble 


(3.3D 


(3.32) 
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rate  by  the  empirical  expression 


0.145037  -  0.56877  x  IO"*t  rtv/s«c  (14,868  <  t<  15^420  stc  UT)  (3. 33) 

where  t  is  measured  in  seconds  after  UT  August  27,  1958. 

The  tumble  phase  angle  Y  (Sec.  3.1,  Fig.  3.1)  was  determined  by 

y*a  rcos(A^),  (3.,4) 

where  is  a  unit  vector  in  Vernal  Equinox  coordinates  (Pig.  3.1). 

The  time  tQ  when  BQ  (magnetic  flux  density  unit  vector  at  tQ,  Pig.  3.1) 
is  in  the  Xg—Yg  plane,  is  then  the  initial  time  at  which  the 
Satellite  Inertial  Coordinate  System  was  determined,  and  t  in  Eq.  (3.11) 
is  expressed  in  seconds  after  t^. 

Naumann's  (1961 )  method  of  determining  the  spin  or  roll  angular 
velocity,  or  rate,  was  used  for  Pass  454  and  attempted  for  Pass 

414.  Equation  (3.11)  Indicates  that  for  times  when  cos  (  ^  &)»  °* 

the  count-rate  of  the  directional  scintillation  counter  will  not  show 
any  variation  due  to  tumble  rotation  and  will  be  modulated  only  by  the 
spin  or  roll  of  the  satellite.  By  the  determination  of  the  period  of 
such  modulations  Naumann  was  able  to  measure  the  time  history  of  the 
roll  period  as  shown  by  the  circle  points  of  Pig.  (3.3).  The  interval 
August  24-29,  1957  exhibited  the  anomalous  Increase  -  decrease  in  roll 
period  (the  Increase  means  a  decreasing  spin  rate,  the  decrease  which 
followed  means  an  Increasing  spin  rate. 

In  the  present  work,  all  attempts  to  use  the  roll  rates  determined 
by  Neumann  for  Pass  414  (August  27,  1958)  failed  to  yield  any  coherent 
picture  of  the  variation  of  the  Detector  A  count-rate  with  the  angle  0 
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(between  the  counter  axis  and  a  plane  perpendicular  to  the  geomagnetic 
field  direction).  Although  all  values  of  6  (in  2°  steps)  were  used 
to  produce  plots  of  d  as  a  function  of  time,  using  Eq.  (3.11),  for  roll 
periods  from  30-60  secs,  no  plot  of  Q  vs  time  could  be  produced  which 
matched  the  variations  in  counting  rate  for  more  than  a  minute. 

A  more  recent  determination  of  the  roll  period  time  history  was 
achieved  by  Bock  (1964)  which  resulted  in  the  non-anomalous  dashed  line 
of  Pig.  (3.5).  Use  of  the  monotonically  increasing  roll  periods  (grad¬ 
ual  slowing  down  of  satellite  spin)  of  Bock  resulted  in  a  coherent  plot 
of  the  Detector  A  count -rate  vs  6  for  Pass  4l4  (see  also  Sec.  3.4) 

Purther  evidence  for  the  correctness  of  Bock's  non-anomalous  roll 
periods  was  developed  in  the  present  work  through  Pourier  series  analy¬ 
ses  of  the  Detector  A  count-rates  for  Passes  414,  427  end  453,  as  is 
detailed  in  Appendix  1,  Power  Spectrum.  This  analysis  Is  concerned  with 
beat-frequency  periodicities  appearing  In  the  count-rate  data.  The 
squares  used  in  Pig.  3.5  show  the  results  of  this  kind  of  analysis  for 
Passes  414,  427  snd  453  and  suggest  strongly  that  the  anomalous  roll 
periods  plotted  by  Neumann  for  the  Interval  August  24-29,  1937  are  due 
to  measurement  of  beat  frequencies  rather  than  single  modulation  fre¬ 
quencies  In  the  Detector  A  count-rates. 

Por  Pass  4 14  Bock  (1964)  evaluated  the  true  roll  rate  by  the  empiri¬ 
cal  expression 

+  -  0.093858  *  0.3103975  x  I0*t  ♦  0.209997%  x  O*  t*  rtv /—c  (3.33) 

where  t  Is  measured  In  seconds  after  the  acquisition  time  for  the  pass 
or  14,868  sec  UT. 


37 


The  roll  phase  angle  6  is  the  most  difficult  parameter  to  deter¬ 
mine  in  Eq.  (3.11).  Naumann's  (196])  method  does  not  work  when  roll 
and  tumble  rates  are  nearly  the  same.  Hence,  6  was  determined  in 
the  present  work  (Program  3-3)  by  a  brute-force  variational  procedure 
termed  "closure-search"  as  follows. 

If  all  the  variables  of  Eq.  (3.11)  except  5  are  determined  0  can 
be  calculated  for  any  time  of  interest  by  performing  a  variation  on  6. 
The  procedure  adopted  was  to  determine  the  time  of  an  r-f  null  shortly 
after  ground  station  acquisition  of  the  satellite.  This  time  was  used 
as  reference  time  tQ  from  which  t  in  Eq.  (3.11)  was  measured.  The 
satellite  inertial  axes  (Sec.  3.1.  Pig.  3.1)  were  calculated  using 
Eqs.  (3.12)  and  (3.13)  and  the  value  of  B  at  tQ.  A  was  computed  using 
Eq.  (3.13)  and  phase  angle  y  as  explained  in  Sec.  3.1. 
Polynomials  (3rd  or  less  order)  which  fitted  the  counting  rate 
for  the  short  (.06  min,  approximately)  intervals  between  signal  fade- 
outs  were  used  to  evaluate  the  Detector  A  count -rates  for  all  desired 
times.  0  was  then  calculated  for  each  of  these  count-rate  determina¬ 
tions  for  all  values  of  6  in  the  range  0  -  360°.  For  each  delta, 
Cobs  aS  a  function  of  X,  j £-  9 1  was  considered  by  the  com¬ 

puter  throughout  the  time  Interval  of  the  polynomial  in  question.  If 
the  polynomial  fitted  data  yielding  the  same  count  rate  at  two  differ¬ 
ent  times  (for  example,  polynomial  e  Pig.  2.10),  ten  time-pairs  of 
equal  count-rates  were  determined  from  the  polynomial  and  pairs  of  X 
corresponding  to  the  time  pairs  calculated  using  Eq.  (3.11).  If  more 
than  70  percent  of  the  X  P*lr*  were  within  3°  of  being  equal  a  graph 
of  the  function  C  b  vs  \  was  plotted.  If  the  polynomial  fitted 
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data  in  which  most  of  the  counting  rate  values  were  unique  (for 

example,  polynomial  c  Pig.  2. loathe  slope  and  intercept  of  a 

straight  line  fitted  to  C  .  vs  X  was  listed  for  each  6  .  This 

obs 

procedure  was  carried  out  for  98  polynomials  for  Pass  414  and  one 
degree  steps  (0-360°)  for  5  ,  including  both  the  natural  and  the 
ARGUS  trapping  regions. 

The  resulting  plots  and  listing  of  slopes  then  were  investigated 

to  determine  a  single  value  of  6  which  satisfied  the  criterion  that 

the  graph  of  Cobs  vs  X  would  be  single  valued  for  each  of  the  98 

sections  of  the  Detector  A  count-rate  intervals  considered.  If  any 

of  the  parameters  of  i£q.  (3.11)  were  slightly  incorrect  the  plots  of 

Cobs  vs  ^  were  *’oun<1  t0  be  multivalued  so  that  Cobg  vs  X  points 

formed  an  open  loop.  The  procedure  as  outlined  therefore  was  termed 

•  • 

a  "closure-search"  and  the  values  of  tumble  rate,  roll  rate  and 

bl,  b2,  b3,  were  considered  to  be  correct  when  a  single  value  of  6 

satisfied  the  closure  criterion  for  all  polynomials  used. 

Por  Pass  414  an  r-f  null  was  observed  at  248. WJ  UT  August  27,  1938. 

Y  was  calculated  to  be  35i*6and  &  was  evaluated  by  closure- 

search  as  90°.  This  value  of  6  satisfied  the  closure  criteria  for 

all  98  sections  of  the  Detector  A  count -rate  data  except  those  in  which 

the  background  counting  rate  varied  rapidly  during  a  time  period  of 

about  one-half  the  tumble  period. 

An  additional  check  on  the  determination  of  parameters  Is  demon- 

4 

st rated  by  Pigs.  2.9  and  2.10  wnereln  Cos  6  vs  time  (Upper)  is  compared 

4 

with  CQ^8  vs  time  (Lower)  for  the  scintillation  counter.  The  Cos  6, 
suggested  by  Bock,  is  an  empirical  analytic  fit  to  the  Coto#(9)  vs  9  curves 
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of  Sec.  2.3  of  the  present  report.  It  is  observed  in  Figs.  2.9,  2.10 
4 

that  cos  6  vs  time  curves  represent  well  the  C  ,  vs  time  curves  for 

obs 

both  the  natural  and  the  ARGUS  trapping  regions  except  for  the  time 
interval  during  which  there  is  sudden  increase  in  omnidirectional 
background  radiation  (curve  B)  in  the  ARGUS  region.  This  background 
radiation  effect  is  discussed  in  Chap.  5* 


l 
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3.4  Observed  directional  count-rate  for  ARGUS  trapping,  C  vs  6. 

OD5 

Results  and  conclusions. 

As  noted  and  explained  in  preceding  sections  of  the  present  chapter, 
the  count-rate  of  the  directional  detector  (Detector  A,  Channel  2)  was 
determined  as  a  function  of  the  angle  6  between  the  scintillation 
counter  axis  and  the  plane  perpendicular  to  the  geomagnetic  field  for 
both  the  natural  and  the  ARGUS  trapped  radiation. 

Figures  2.9  (natural  trapping,  Pass  4l4)  and  2.10,  2.11  (ARGUS  trap¬ 
ping,  Pass  .414)  illustrate  the  results  of  the  analyses  of  directional 

4 

count-rate  (CQl}S)  and  angle  0  (in  form  of  cos  0)  vs  real  time  (UT). 

Figures  3.4,  3.5  and  3.6  illustrate  the  results  of  the  determination 
(Sec.  3.3,  closure-search  procedure)  of  observed  directional  count-rate 
C'obs  vs  0  for  several  of  the  modulations  designated  a,  b,  c  ....  q  in 
Figs.  2.10  and  2.11. 

Although  the  "disc -like"  character  of  the  trapped  radiation  is  gener¬ 
ally  evidenced  by  the  maxima  of  CQb3  at  0  =  0°  and  decreased  values  at 
larger  angles,  the  plots  of  Figs.  3.4,  3.5  and  3*6  show  considerable 
variation  in  the  angular  distributions  of  directional  count-rates 
(CQb8  vs  0)  for  ARGUS  trapping.  These  variations  are  considered  to  be 
due  mainly  to  1.  scatter  in  the  data  Itself,  and  2,  effects  of  the  omni¬ 
directional  penetrating  radiation  (curves  B,  Pigs.  2.10  and  2.11)  upon 
the  counting-rate  of  directional  Detector  A. 
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Item  2  just  noted  is  considered  to  be  one  of  the  special  conclusions 
of  the  present  work.  However,  since  it  has  been  possible  in  the  time 
available  to  analyse  but  a  small  percent  of  the  available  Explorer  IV 
directional  count-rate  data  by  the  closure-search  procedure  (ijec.  3.3), 
any  conclusions  stated  here  must  be  considered  as  tentative.  Keeping 
this  reservation  in  mind,  there  is  suggestive  and  fairly  convincing 
evidence  nevertheless  that  the  shape  of  the  angular  distribution  of 
count-rate  (CQbs  vs  6,  Pigs.  3.4,  3. 5  and  3.6)  is  affected  by  the  omni¬ 
directional  penetrating  radiation  (curves  B,  Pigs.  2.10  and  2.11);  as 
follows : 

(1)  Plot  b,  c,  d  (Pig.  3.4,  upper  left)  is  a  combined  plot  of  modu¬ 
lations  b,  c,  d  of  Pig.  2.10,  Lower.  Thus  all  3  modulations  can  be 
represented  in  a  single  plot,  suggesting  that  the  low  intensity  and 
constant  background  penetrating  radiation  (curve  B;  Fig.  2.10,  Lower) 
is  not  a  disturbing  factor  in  the  Co^8  vs  6  analysis. 

(2)  Plot  f  (Pig.  3.4,  upper  right)  is  an  example  of  the  apparent 
effect  of  the  rapidly  increasing  background  (curve  B,  Fig.  2.10 
Lower),  it  being  suggested  that  the  set  of  points  indicating  abnor- 
aally  large  count-rates  (near  6  »  J>0°-60°)  Include  some  background 
radiation.  It  appears  that  the  CQ^s  vs  6  count-rate  data  should  be 
corrected  for  changing  background  as  the  satellite  moves  through  the 
AAQUS  belt  for  regions  such  as  f  (Pigs.  3.4,  2.10). 

(3)  Plot  g  (Pig.  3.4,  lower  left)  exhibits  a  broader  angular  distri¬ 
bution  of  the  directional  count-rate.  At  the  same  time  curve  B,  Pig.  2.10 
Lower,  Involves  an  increasing  background  penetrating  radiation  which  may 
have  contributed  to  the  directional  count-rate  intensity  level  to  cause 
the  broadening  of  curve  g.  It  is  difficult  to  make  a  quantitative  esti- 
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mate  of  the  precise  contribution  of  the  omnidirectional  G-M  count-rates 
to  the  directional  scintillation  detector  count-rates  since  the  calibra¬ 
tions  of  the  detectors  involved  did  not  include  intercalibration  of  the 
G-M  and  scintillation  counters  for  the  same  penetrating  radiations. 

Also,  any  penetrating  component  of  the  directional  radiation  might  be 
more  readily  detected  by  the  scintillation  counter  at  larger  8  values 
since  there  is  generally  more  shielding  in  directions  away  from  the 
counter  axis. 

(4)  Plot  h  (Pig.  3.4,  lower  right)  has  broader  angular  distribution 

than  plot  g,  the  penetrating  background  radiation  curve  B,  Pig.  2.10, 

Lower,  is  at  the  same  time  somewhat  more  Intense  than  for  plot  g  and 

so  the  effect  of  the  penetrating  background  radiation  is  quite  similar 

to  that  suggested  for  plot  g.  Further  the  count-rate  of  the  directional 

4 

detector  is  approaching  saturation,  the  rate  being  in  excess  of  10  per 
sec,  a  situation  which  would  produce  relative  broadening  of  the  angular 
distribution  by  the  flattening  process. 

(5)  Plot  J  (Pig.  3-5,  upper  left)  corresponds  to  a  decreasing  back¬ 
ground  radiation  (curve  B,  Pig.  2.11)  suggesting  that  the  narrowness  of 
the  angular  distribution  of  plot  J  is  due  to  decreasing  directional 
count-rate  at  larger  6  values. 

(6)  Plots  k,  1,  a,  n,  o  (Pig.  3*5  and  3.6)  show  variation  in  the 
angular  distribution  of  the  directional  count-rates  similar  to  the  plots 
of  Pigs.  3.4  and  3*5  already  discussed  and  these  variations  suggest  again 
the  effect  of  penetrating  background  radiation  (curve  B,  Pig.  2.11).  Plots 
a  and  o  suggest  that  a  burst  of  penetrating  radiation  was  encountered 
toward  the  end  of  the  tine  periods  for  the  two  plots  since  the  two  groups 
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of  points  showing  maximum  count-rate  for  ranges  of  6  of  about  0°- 
30°  or  40°  occur  last  in  the  time  sequence.  This,  penetrating  radia¬ 
tion  effect  may  be  present  for  the  entire  sequence  of  plots  1  through 
p.  The  effects  are  not  always  clear  and  unmistakable  but  they  do 
appear  to  be  present. 

(7)  Finally  plots  p  and  q  (Fig.  3.6  lower)  seem  to  show  the  effect 
of  the  sudden  decrease  in  curve  B,  Fig.  2.11,  occurring  at  about  the 
same  time  as  the  sharp  decrease  in  count-rate  for  plot  p  (near  6  =  5°) 
and  the  continued  less-sharp  decrease  for  plot  q.  These  sudden  de¬ 
creases  in  count-rates  are  similar  effects  to  the  corresponding  in¬ 
crease  shown  by  plot  f. 
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Table  3.1.  Portion  of  1958  Epsilon  (Explorer  IV,  Pass  4l4)  Orbital  Data  Series 
(Smithsonian  "Astrophysical  Observatory,  1959) 


Time, 

UT 

(hr-min) 

East 

latitude 

(°) 

Geo¬ 

graphic 

Latitude 

(°) 

Height 

(km) 

Geomagnetic  Flux 
Density  Components 
(Gauss) 

Satellite  Velocity  Components 
(kiVhr) 

bn 

be 

bd 

VN 

VE 

VD 

04 

05 

268.34 

-06.85 

1953.4 

.136 

0.019 

0.019 

014996.1 

16001  .8 

01  745.6 

Ok 

06 

269.84 

-04.76 

1923.1 

.138 

0.018 

0.031 

01  5089.3 

16035.6 

01  874.0 

Ok 

°z 

271.34 

-02.64 

1890.7 

.140 

0.016 

0.043 

01  5173.0 

16089.5 

01999.1 

04 

08 

272.85 

-00.51 

1856.3 

.141 

0.015 

0.056 

01  5246.8 

16163.6 

02120.7 

04 

09 

274.38 

001 .64 

1820.0 

.142 

0.013 

0.069 

015310.4 

16258.5 

02238.4 

04 

10 

275.92 

003.81 

1781.8 

.144 

0.01 1 

0.082 

01  5363.2 

16374.4 

02351  .8 

04 

11 

277.49 

005.99 

1741.9 

.144 

0.009 

0.096 

01  5404.1 

16512.2 

02460.6 

04 

12 

279.09 

008.20 

1700.2 

.145 

0.007 

0.110 

01  5432.1 

16672.5 

02564.4 

04 

13 

28O.73 

010.41 

1656.9 

.145 

0.005 

0.125 

01  5445.7 

16856.1 

02662.8 

04 

14 

282.43 

Q12.64 

1612.0 

.146 

0.002 

0.139 

01  5443.1 

17064.0 

02755.5 

04 

15 

284.17 

01 4.88 

1565.8 

.145 

-.001 

0.154 

01  5422.0 

17297.2 

02842.0 

04 

16 

285.99 

017.13 

1518.2 

.145 

-.004 

0.170 

015379.7 

17  557.1 

02921 .8 

04 

17 

287.88 

019.38 

1469.4 

.144 

-.008 

0.185 

015313.2 

17844.8 

02994.4 

04 

18 

289.86 

021 .64 

1419.4 

.142 

-.012 

0.200 

01  5218.3 

18162.0 

03059.5 

04 

19 

291.94 

023.89 

1368.6 

.140 

-.016 

0.215 

01  5090.7 

18510.3 

03116.6 

04 

20 

294.13 

026.13 

1316.9 

.138 

-.020 

0.229 

014924.8 

18891 .0 

03165.0 

04 

21 

296.45 

028.36 

1264.4 

.136 

-.024 

0.243 

014714.1 

19306.2 

03204.4 

04 

22 

298.92 

030.57 

1211.5 

.134 

-.029 

0.256 

014450.9 

19757.1 

03234.2 

04 

23 

301.55 

032.75 

1158.1 

.131 

-.033 

0.268 

014126.4 

20245.0 

03253.9 

04 

24 

304.38 

034.89 

1104.5 

.129 

-.038 

0.279 

01 3730.3 

20770.5 

03262.9 

04 

25 

307.40 

036.98 

1050.9 

.127 

-.042 

0.289 

01 3250.5 

21333.7 

03260.8 

04 

26 

3!  0.67 

039.00 

0997.4 

.126 

-.046 

0.298 

01 2673.6 

21933.0 

03247.1 

04 

2Z 

314.18 

040.93 

0944.2 

.125 

-.049 

0.305 

011985.2 

22564.5 

03221  .4 

04 

28 

317.97 

042.76 

0891.5 

.125 

-.052 

0.312 

011169.4 

23222.7 

031 83.1 

04 

29 

322.07 

044.46 

0839.5 

.125 

HEEJ 

0.317 

01 0210.9 

23897.3 

031 31.9 

04 

30 

326.47 

046.01 

0788.3 

.127 

BGE3 

0.321 

009095.9 

24574.4 

03067.4 

04 

31 

331.20 

047.37 

.130 

-.053 

0.325 

007814.7 

25234.3 

02989.4 

04 

32 

336.24 

048.52 

0689.7 

.133 

-.050 

0.328 

006362.7 

25853.8 

02897.7 

04 

33 

341.58 

049.43 

0642.6 

.138 

-.046 

0.330 

004745.8 

26404.7 

02792.1 

04 

34 

347.17 

050.07 

0597.3 

.143 

-.040 

0.333 

002982.2 

26858.0 

02672.7 
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Table  5.2.  Orbital  data  of  Table  3.1.  with  magnetic  flux  density  components 
in  Vernal  Equinox  Coordinate  System 


Time, 

UT 

(min) 

East 

Longi¬ 

tude 

(°) 

Geo¬ 

centric 

Lati¬ 

tude 

(0) 

Distance 
Satel¬ 
lite  to 
Geo¬ 
center 
(Earth 
radii) 

Components  of  Mag¬ 
netic  Flux  Density 
Unit  Vector  in  VE 
Coordinate  System 

Satellite  Velocity  Components 
(km/hr) 

BX 

H 

m 

VN 

mm 

245. 

246. 

247. 

248. 

249. 

250. 

251. 

252. 

253. 

254. 

255. 

256. 

259. * 

260. 
261. 
262. 

263. 

264. 

ill: 

18: 

269. 

270. 

271. 

272. 

St 

268.34 

269.84 

271 .34 

272.85 

274.38 

275. 92 
277.49 
279.09 
280.73 
282.43 

284.17 
285.99 
287.88 

289.86 
291 .94 
294.13 
296.45 

298.92 
301.55 

304.38 
307.40 
310.67 

314.18 
317.97 
322.07 
326.47 
331.20 
336.24 
341.58 
347.17 

-6.80 

-4.72 

-2.62 

-.50 

1  .62 
3.78 

$•95 

8.14 

10.34 

12.55 

14.78 
17.02 

19.25 
21 .50 
23.74 
25.97 
28.19 
30.40 
32.57 
34.70 

36.79 
38.81 
40.73 

42.56 

44.26 
45.81 
47.17 
48.32 
49.23 
49.87 

1  .3062 
1  .3014 
1 .2964 
1  .2910 
1  .2853 
1 .2793 
1  .2730 
1.2664 
1.2596 
1  .2525 
1  .2452 
1  .2377 
1 .2300 
1 .2220 
1 .2140 
1  .2058 
1.1974 
1.1890 
1.1805 
1.1720 
1.1635 

\:\m 
1.1382 
1.1299 
1.1218 
1.1139 
1 .1062 
1 .0988 
1.0916 

2 

W  4 

•  f A 

f  ^ 

.093 

.185 

.262 

.339 

.399 

.445 

.486 

.512 

.531 

.532 

.527 

.512 

.486 

.453 

.413 

.367 

.316 

.257 

.197 

.129 

.062 

-.008 

-.°7j 

-‘.220 

-.291] 

-.359 

-.423 

-.W 

-.546 

.990 

.982 

.962 

.928 

.883 

.832 

.769 

.702 

.627 

.556 

.477 

.397 

.319 

.239 

.161 

.088 

.018 

-.048 

-.112 

-.170 

-.224 

-.271 

-.313 

-.349 

-.382 

-.405 

-.423 

-.438 

-.443 

-.445 

14996.1 

15089.3 
15173.0 
15246.8 

15310.4 

1  5363.2 

1  5404.1 

1  5432.1 

15445.7 

1  5443.1 

1  5422.0 

1  5379.7 

15313.2 

1  5218.3 

15090.7 

14924.8 

14714.1 

14450.9 

141 26.4 

13730.3 

13250.5 

12673.6 

11985.2 

11169.4 

10210.9 
9095.9 

7814.7 

6362.7 

4745.8 
2982.2 

16001  .8 
16035.6 

16089.5 

16163.6 
162  58.5 

16374.4 
16512.2 

16672.5 
16856.1 
17064.0 

17297.2 

ipi 

1 81 62.0 

18510.3 
18891.0 
19306.2 
19757.1 
20245.0 

20770.5 

21333.7 
21933.0 

22564.5 

23222.7 

23897.3 

24574.4 
25234.3 

25853.8 
26404.7 
26858.0 

1745.6 
1874.0 

1999.1 

2120.7 

2238.4 
2351  .8 
2460.6 

2564.4 

2662.8 

2755.5 
2842.0 
2921  .8 

2994.4 

3059.5 

3116.6 
3165.0 

3204.4 

3234.2 

3253.9 

3262.9 
3260.8 

3247.1 
3221  .4 

3183.1 

31  31.9 

3067.4 

2989.4 

2897.7 

2792.1 

2672.7 
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Pig.  3.1  Geometries  employed  in  determination  of  count-rates  of 
Explorer  IV  directional  scintillation  counter  in  terms  of 
orientation  of  the  satellite  relative  to  geomagnetic  field 
direction.  All  symbols  are  explained  in  the  text. 
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SATELLITE  INERTIAL  COORDINATE  SYSTEM  SATELLITE  FIXED  COORDINATE  SYSTEM 


Pig.  3. 


2  Sketch  showing  various  geometric  and  vector  quantities 
and  their  interrelations.  A,  geodetic  latitude; 

Zi  ,  geocentric  latitude. 
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Pig.  3.3  Roll  or  spin  period  time  history  of  Explorer  IV  Satellite 
showing  anomalous  increase-decrease  (circle  points)  between 
August  24  and  29,  1958  -  calculated  by  Naumann  (1961). 

Dashed-line  (X  points)  represent  a  non-anomalous  roll 
period  history  as  determined  by  Bock  (1964)  by  an  analog  computer 
method. 

Squares  locate  beat -frequency  periods  measured  in  the 
Detector  A  count-rate  data  by  the  authors  of  the  present  report 
and  suggest  that  Naumann* s  anomalous  roll  periods  are  due  to 
measurement  of  beat  frequency  periods  rather  than  single  (correct) 
roll  periods. 
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Pigs.  3.4,  3.5,  3.6  Observed  directional  count -rate  CQba  vs  9  during 
penetration  of  ARGUS  (Event  I)  shell.  Pass  4l4. 

The  lower  case  letters  in  parentheses  in  the  upper  left  hand 
corners  of  the  individual  plots  correspond  to  the  same  lettered 
count-rate  modulations  shown  by  Pigs.  2.10  and  2.11.  The  twelve 
separate  plots  were  obtained  by  individual  closure-search  proce¬ 
dures  as  explained  in  Sec.  3*3. 

Abscissa  values  for  6  were  calculated  by  Eq.  (3.11)  using 
observed  times  t  and  values  of  the  other  parameters  as  discussed 
in  the  text.  The  same  time  references  (248.941  UT),  and  phase 
angles  (  y  *  351.6°,  S  «  50°)  were  employed  for  all  the 
plots  of  Pigs.  3.4,  3.5,  3.6. 

The  various  steps  in  the  reduction  and  analysis  of  data  followed 
to  obtain  the  data  for  CQbg  vs  $  (Pigs.  3.4,  3.5.  3.6)  are 

(1)  CQb8  vs  time  is  obtained  as  in  Pigs.  2.10,  2.11; 

(2)  Polynomial-fitting  to  the  plots  of  (1)  is  accomplished; 

(3)  y-  and  S  are  evaluated  from  r-f  null  data  and  the 
closure-search  procedure  respectively; 

(4)  b. ,  b  and  b  are  calculated  from  B  ,  B  ,  B„  In  the  Vernal 

1.  c  j  x  y  z 

Equinox  Coordinate  System  and  Neumann's  estimate  of  angular 
A 

momentum  L  direction. 

/V  #  • 

(5)  The  magnitude  of  L,  ^  and  ^  are  obtained  using 
J.  Bock's  decay  coefficients; 

(6)  6  Is  calculated  by  Eq.  (3.11). 
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(«e-l)  Wsec 


D»S)Sq<b 


Chapter  4 


Determination  of  the  Directional  Flux  Density  T(e')  from  Observed 
Directional  Counting  Rate  C(®) 


4.1  Conversion  of  the  integral  equation  to  matrix  form 

The  directional  flux  density  is  related  to  the  directional  counting 
rate  by  the  integral  equation 


C«- 


v  v 


SMT(®’)ari 


where 


O) 


-2  -1 

is  the  true  counting  rate  (cm  -sec  )  expressed  as  a 


function  of  the  angle  between  the  geomagnetic  flux  density  direction  B 
(Pig.  4..1)  and  the  direction  of  the  scintillation  counter  (Detector  A) 

axis  P,  S(*0  «  S  f)  is  the  Detector  A  response 

A 

function  (  e(j  angle  between  the  direction  of  the  counter  axis  P  and 
the  unit  vector  which  defines  the  orientation  of  the  element  of  solid 
angle Pig.  4.1)  and  is  the  directional  flux  density 

(partlcles/cm2sec-ster).  The  observed  counting  rate  ( ®)  18 

corrected  for  detector  deadt  lines  to  yield  the  true  counting  rate 
c<®)  by  the  equation: 

C*.« 


C(®)  *  £ 


A  l-nC*.(«) 


(4.1) 


(4.2) 
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A  is  the  effective  area- of  Detector  A  (0.005  cm^;  Baicy,  1 96 2 ) ;  r, 
is  the  detector  deadtime  which  has  been  estimated  as  100  microseconds 
(Chapter  2).  Prom  Pig.  4.1  K  is  related  to  8^  <f>  *  and  0  by 


Cos»<  *  Cos0,5i'neSt'n^,-i-  Sin e‘  Cos 6 


The  Detector  A  response  function 


SW 


is  shown  by  Pig.  4.2  (Baicy, 


(4.3) 


1962).  A  good  analytic  fit  to  the  curve  of  Pig.  4.2  is  (Weber  et  al, 

1962) 

SW  =  e/j>( >  (4.4) 

where  Yj  =  -.00274,  ” *000006151*  and  oC  is  expressed 

in  degrees.  In  the  determination  of  U'(e')  9  integral  Eq.  (4.1) 
must  be  converted  to  matrix  form  for  computer  purposes.  To  accomplish 
the  evaluation  of  the  integral  of  Eq.  (4.1)  in  selected  intervals  of 
6j  j(eO  is  assumed  constant  in  the  interval .  por  such 
step-by-step  Integrals  the  limits  of  Integration  of  0*  and 

are  determined  from  an  analysis  of  the  possible  orientations  of  the 

A 

counter  axis  with  respect  to  the  magnetic  flux  density  B  and  the  parti¬ 
cle  flux. 

Also  cylindrical  symmetry  of  jj(8^  relative  to  the  cylinder  axis  in 

A 

the  direction  of  the  magnetic  flux  density  B  was  assumed  so  that 

J(«')  ■  T(-«')  (  where  J  (Q^ 

is  Independent  of 

)•  Prom  the  response  function  Pig.  4.2  it  is  observed  that 

a  cone  of  half -angle  30°  approximately  and  with  its  axis  along  the  counter¬ 
axis  approximately  defines  the  region  in  which  particles  will  be  counted. 
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However,  for  the  purposes  of  analysis,  a  hemisphere  instead  of  a  cone 
is  considered  for  such  a  geometrical  construct  does  not  deviate  appre¬ 
ciably  from  the  actual  physical  situation  because  of  the  negligible 
value  of  S(*0  for  >  30°.  Prom  Pig.  4.1  it  can  be  seen  that 

the  limits  of  integration  of  p  depends  in  general  on  0*  and  0  ; 

takes  on  all  values  from  0  to  2TT  for  some  values  of  0*  and 
goes  from  0  +  ^  '  to  2 TT  —  Alp*  for  other  values ;  where 

A</>'  >  O  and  depend  on  0*  and  0, 

Prom  geometrical  considerations,  it  can  be  shown  that  A<p  m 
Arcsin  (  tan.  0  *  tan,  O)  •  In  Eq.  (4.5)  O  ^  0  **  and 

4  0  ^  •  Therefore  ( IflU  0 *  IflTU 0 )  can 

be  greater  than  1  for  certain  values  of  0*  and  0  .  Since 

has  to  be  real,  Eq.  (4.5)  is  restricted  to  those  values  of 
0  *  and  0  for  which  T5tV  0  fan,  0^1 
the  value  of  for  (  fan  0*  tun,  0  )  >  1 

It  is  shown  readily  that  Eq.  (4.5)  is  valid  in  the  fonn 

Af'  a  /4rctm(T&7teT<m«&,(®+8)sf,  (4.5) 

rime  ©'+  0  =  0. 

The  analysis  above  shows  that  the  limits  of  integration  for  P 
are  from  0  to  2  IT  when  0*  and  0  are  such  that  e‘+e>J  ; 
and  when  8*+  0  ^  ^  y  $  ranges  from  A<j>‘ 

to  TT-  Af  for  0*  varying  from  -  (f-e)  to  OJ 

^  ranges  from  to  TT  +  4^*  for  0 

Thus  Eq.  (4.1)  can  be  written  in 


.  Also, 
is  required. 


varying  from  0*  to  ZE  —  0 

a 

the  explicit  form 
,0 


C{9) 


l 


¥-0  ir+Ap'  it  tir 

r*  r  /» X  r 


V  ®*o  it*  tff-*  Jo 


(4.6) 
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To  approximate  the  integral  of  Eq.  (4.6)  by  a  matrix  equation  it  is 
assumed  that  1(f)  varies  slightly  only  when  0  is  varied 
by  ±  2.5°.  Thus  1(f)  is  required  only  for  values  of  0*  given 
by  (  A.  “  4  *  I)  •••  1  8  j  and  *  <5  •  The 

value  of  J(e‘)  anywhere  inside  the  interval  of  width  5°  is  given 
by  the  linear  approximation 

<•-«> 


LH 


(4.7) 


where 


9[<e'<  e!  . 


A  similar  type  of  reduction 
is  carried  out  for  the  C(e)  function.  Thus 

Cjl  3  C  (0>)  }  where  0^  -  40  and  ^0*5* 

Equation  (4.6)  can  now  be  written  as 


M  C  r 

r‘C  ^ 


I  +  J  J 

V*  -Ap1  J 


11  fV*1  r 

+  Z  T(v) 


2TT 

swan' 


(4.8) 


where  and  the  value  of  the  Integer  M  Is 

I  *nr 

to  be  chosen  such  that  0  +  9:  &  t  Now,  noting  that 

4  J  • 

0*)  *  0"  (“0')  and  also  assuming  that  T(0') 

to  be  linear  In  a  sub-interval  of  width  5*.  the  Integral  Eq.  (4.8) 
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reduces  to 


M  n>~  & 


cta)«Q  =na; 


M 


k*1 


rl  .V,  ^ 

^>Tr+A^»i 

5.Wan'+ 
&&  *  rM+t 


-rr-  Afyt 

SMlri  + 


*V«  ^7^ 


^(<)JA  5 


where  djfi  =  CosG'de'd^1  , 

A *  /Arcsia  ( tan  0,1.  tan©;)  , 

5X to*  S(^,e\r). 


(4.9) 


It  can  be  shown  that  the  value  of  M  depends  on  the  index  1  and  is 
given  by  18-  i.  The  number  18  is  the  number  of  sub-intervals  of 
width  5#  in  the  interval  0*<  ©'  30*  •  Equation  (4.9)  can 

be  written  in  the  matrix  form 

18 

Ql  *  S* j  Dj  y  t,J  s  1^....  18  j  (4.10) 


where 


and  S^j 


or- 


J(ej  ) ,  , 

c  rv^? 

n  si-)^ .  for  i+j  ^  16  . 
7  O 
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Thus  the  reduction  of  the  integral  form  of  Eq.  (4.1)  to  the  matrix 

form  of  Eq.  (4.10)  is  completed  when  all  the  324  elements  of  the 

S-matrix  are  determined.  These  matrix  elements  have  been  computed 

using  the  techniques  of  numerical  integration.  If  T(o')  is  con- 

- 1 

sidered  independent  of  “  *  the  counting  rate  would  be  Independent 

l8 

of  direction  of  the  counter  axis  and  hence  the  sum  2  S  £ { 

J*l  J 

should  be  independent  of  the  index  i.  The  computed  values  of  the 
S-matrix  elements  have  been  shown  to  satisfy  this  condition  with  con¬ 
siderable  accuracy  (  per  cent  differences  are  of  the  order  of  .025 %  )• 

4.2  Inversion  of  the  matrix  equation.  Equation  (4.10) 

The  determination  of  involves  solving  18  simultaneous  equa¬ 
tions,  Eqs.  (4.10).  For  a  given  set  of  9  a  unique  solution 

exists  if  the  det(S)/  zero.  In  the  present  case,  the  S-matrix  ele¬ 
ments  have  the  physical  interpretation  of  an  approximate  measure  of 
the  resolving  power  of  the  directional  detector.  The  det(S)«zero  has 
the  physical  meaning  that  the  scintillation  counter  would  have  zero 
directional  sensitivity.  From  the  shape  of  the  response  function,  the 
determinant  of  S  should  be  non-zero.  However,  it  has  been  found  that 
det(S)fl<10~^.  However,  it  has  been  possible  to  Invert  the  S-matrix 

with  sufficient  accuracy  using  an  IBM  1620  computer  so  that  numbers 
-99 

less  than  10  '  only  are  treated  as  zero.  Mathematically,  the  problem 

is  straightforward  once  the  inverted  S-matrix  elements  are  computed 
since  the  3^  can  then  be  determined  by  the  matrix  equation 

m  2  (S*1). .  Cx  >  i»l,a,...IS.  (4.111 


56 


The  accuracy  of  the  matrix  inversion  has  been  tested  by  computing  the 
product  * S  S^j, S),j  an(i  comparing  it  with  the  identity  matrix. 

This  comparison  shows  that  for  an  assumed  set  of  3"*  used  to  com¬ 
pute  a  set  of  C:  by  Eq.  (4.10),  the  computed  values  of 
determined  by  Eq.  (4.11)  are  correct  within  l.Oljj.  This  accuracy  will 
be  maintained  if  the  C:  values  are  known  accurately;  ie,  within  an 

V 

error  determined  mainly  by  the  magnitude  of  the  inverted  S-matrix  ele¬ 
ments. 

If  "measured"  values  of  are  derived  from  the  flux  density 

function  axe';  then  the  actual  values  are  obtained  from  the 

_  ii 

basic  Eq.  (4.1).  If  the  measured  values  of  are  denoted  by 
and  related  to  by  the  equation  j  where  the 

actual  value  of  is  not  known,  then  by  Eq.  (4.11) 

ay  =  S  (s')ki  c "  , 

S£(S\[C1+A;]  , 

=  2  (s*')uAi  +  crk 

S  AT|,  +  CT*  }  S  •  (4.12) 

By  Eq.  (4.12)  it  can  be  seen  that  the  error  In  the  computed  value  of  J 
is  proportional  to  the  absolute  value  of  (s'  )*J  for  a  given  set  of 
.  It  has  been  found  that  the  values  of  the  inverted  matrix 
elements  lie  between  -3000  and  +5000  so  that  an  error  of  lji  in  the  count¬ 
ing  rate  would  give  rise  to  unreliable  J-values.  Although  S*1  has  been 
computed  with  sufficient  accuracy,  owing  to  the  smallness  of  det(S), 
small  errors  in  Cl  can  produce  physically  meaningless  Cjj  ;  this  has 
been  tested  by  using  C  values  with  assumed  values  of  /^^and  it  has 
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been  found  that  the  measured  C-Value:5  are  n°t  accurate  enough  to 
produce  physically  meaningful  J*-  by  Eq.  (4.11).  Hence  it  has 
become  necessary  to  search  for  other  methods  of  solving  ill-conditioned 
matrix  equations.  Recently,  a  new  method  has  been  developed  to  solve 
such  ill-conditioned  equations.  The  method  essentially  involves  the 
construction  of  a  matrix  C  which  satisfies  the  sum  condition  on  S-^ 
and  which  produces  approximately  the  same  result  as  when  very  accu¬ 
rate  data  is  used.  For  this,  it  is  expected  that  the  C-matrix  elements 
will  be  small  in  magnitude  and  tend  to  form  the  matrix  elements  in 
a  certain  limiting  case.  The  computed  C-matrix  elements  do  satisfy 
these  conditions  and  have  been  tested  in  many  cases  of  assumed  T(e') 
functions  and  have  been  found  to  be  satisfactory  for  the  present  problem. 

An  outline  of  the  derivation  of  the  C-matrix  is  given  as  follows: 

For  the  equations 

C,  *  +  S,2  J2. 

C2  =  ^21 

C3  =  Sj,cJJ  +  Sj2^2  (413) 


where  J.  and  are  the  unknowns,  using  only  one  pair  out  of  three  possi¬ 
ble  pairs,  two  unique  values  of  can  be  determined  provided  that  the 
corresponding  coefficient  determinants  are  not  zero.  But  if  a  small 
error  is  introduced  in  the  value  of  such  that  c3  *  =  W 

then  there  are  three  different  pairs  of  J-values 
(OVi  fa**  )  corresponding  to  the  three 

pairs  of  C~values,  (£•!)£*)>  (C,  9  Cj)  J  (Caj  Cj)  .  However,  from 
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the  construct 


3i  \  _  f  )||  (81Y5,'  S^  Sji  V  cA 

jj.y  \ (b  '4 ,  Cb  XftSn.  s»  sj^+h)  (4 14) 

where 


it  can  be  seen  that  the  solution  obtained  from  Eq.  (4.14)  may  be  dis¬ 
tinct  from  all  the  three  pairs  (<T i )  U2.  )>  (  ^1J) 

owing  to  the  fact  that  the  error  in  C3  has  been  averaged  out  in  some 
way  in  producing  the  set  >  31  )•  Also,  it  can  be  seen  that  in 

the  limit  of  £jS  O  *  the  solution  produced  by  Eq.  (4.14)  is  the 
same  as -any  one  of  the  pairs  (a-,1 , 31'  )*  Gr.4>  ?S)  ,  ( T,*,  31»). 
This  method  of  averaging  out  the  error  in  the  data  can  be  generalized 
to  the  case  of  a  square  matrix  by  using  an  intermediate  matrix  which 
has  the  effect  of  averaging  the  errors  In  the  counting  rate  values. 

Thus,  consider 

(Jn-.vj,  % ae'  =  % 

J(e')  -  S  A^^Ce1) * 


M  N 


r 

i 
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so  that  iiq.  (4.10)  can  be  written 


18 


Ci-gs^j;  , 


Next,  define 


SiJ^>  ’ 


Multiplication  of  both  sides  of  I5q.  (4.17)  by  the  transposed  HT  yields 

2  T i,Cx  =  z  E(Ti,T;jAlt-2  B,  A  , 

<«i  Xs i  u.ai  Ay  r  r  as |  yr 


where 


and  multiplication  of  both  sides  of  2q.  (4.18)  by  the  inverse  of  the 
B-matrix,  yields  further 


(4.16) 

(4.17) 

(4.18) 


By  Eq.  (4.15)  substituting 


Ay  from  2q.  (4.19),  yields 


(4.20) 


Equation  (4.20)  becomes,  by  interchanging  the  order  of  summation, 


(4.21) 


Equation  (4.21)  is  of  the  form 


It  can  be  shown  that  in  deriving  the  C -matrix,  the  value  of  N  must  be 
chosen  less  than  the  maximum  value  which  the  i  or  J  Index  can  take 
(here,  it  is  18).  When  N  is  greater  than  18,  the  determinant  of  B 
can  be  shown  to  be  zero  which  renders  the  C-matrix  meaningless.  When 
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N  it  can  be  seen  that  C\j  *  Lj  "or  a11  allov/ed  c^j 

values.  Hence  N  should  be  chosen  less  than  13.  It  is  to  be  noted 
that  for  ^  too  close  to  X  the  initial  approximation  of  a  function 
T(e')  by  a  few  elements  of  an  orthonormal  set  becomes  Inaccurate 
when  the  range  of  0  is  large.  Thus  an  optimisation  of  N  was 
attempted  from  both  theoretical  and  experimental  sides  and  it  was  found 
that  a  good  value  of  N  is  12  for  the  case  of  a  polynomial  of  even 
powers  is  used  for  the  orthonormal  set.  It  has  been  determined  further 
that  the  selection  of  N  should  be  based  on  (1)  the  predicted  error  in 
the  Q- values,  (2)  the  rate  at  which  £  drops  to  its  minimum,  (3) 
the  error  in  recomputing  the  C*  values  from  the  ^".values  computed 
from  Eq.  (4.20),  (4)  the  type  of  orthonormal  set  used  in  computing  the 
matrix.  The  analysis  has  shown  that  in  the  computation  of 
using  Eq.  (4.20)  the  selection  of  h|  based  on  item  (3)  can  be  attained 
by  computing  the  C.  matrix  elements  for  M  #  3 ,  #  £*,  6,7,  ?  and 

using  cosine  functions  as  the  basic  orthonomal  set  of  even  functions. 
Thus,  in  a  more  explicit  form,  Eq.  (4.20)  may  be  written 


j  N  *  If  Zj  •  •  \1  • 


(4.22) 


In  the  final  determination  of 


a  selection  is  made  out  of  the 


set  (  sT  ^  ^  J*  J*  0  9  ,  J*  *  )  based  on  the  accuracy  with 
which  the  input  Q- values  can  be  recomputed  from 


(4.23) 


o2 


It  was  determined  that. in  all  practical  cases,  the  percentage  difference 

T>C  «  aoo(C?-CiVC c?  +  c;)  is  within  5#  which  is  less 

than  the  errror  in  the  £-data. 


4.3  Directional  count-rates  Co^s(0)  and  C  (0)  and  flux  densities  J(0') 
vs  0,  0'  for  Natural  and  ARGUS  trapping.  Results  and  conclusions. 

The  analysis  and  procedure  of  Secs.  4.1  and  4.2  above  were  applied 
to  Passes  4l4  and  454  of  Explorer  IV  over  Huntsville,  Alabama  for  both 
the  Natural  and  ARGUS  trapping  to  attain  observed  count-rates  Co^s(0) 
and  true  count-rates  (observed  count-rates  corrected  for  deadtimes, 
effective  counter  area  and  efficiency)  C(0)  vs  ©(between  scintillation 
counter  axis  and  the  plane  perpendicular  to  the  geomagnetic  field), 
and  directional  flux  densities  J ( 0  1 )  vs  0'.  Figures  4.3,  4.4,  4.5  and 
4.6  along  with  explanatory  captions  are  the  results. 

Any  conclusions  concerning  the  directional  flux  density  distribution 
of  trapped  particles  must  be  considered  as  tentative  since  Passes  4i4 
and  454  only  have  been  analysed  for  J(0')  vs  0'.  The  reason  for  the 
limited  amount  of  data  analysed  to  date  is  the  inherent  nature  of  the 
ARGUS  trapping.  A  great  deal  of  computational  work  is  involved  in  the 
analysis  of  the  data  to  achieve  directional  flux  density  distributions 
for  ARGUS  trapping  as  has  been  noted  in  Chaps.  2,  3.  4.  It  is  necessary  to 
analyse  each  modulation  of  the  observed  count-rate  CQb8  vs  time  (see 
Figs.  2:9*  2.10,  2.11  Lower)  separately  to  attain  a  aeries  of  corres¬ 
ponding  J(0'  )  vs  0*  curves.  In  the  present  work  this  detailed  analysis 
has  been  completed  for  Pass  414  only.  For  the  limited  results  Indicated 
by  Figs.  4.3  -  4.6  the  following  "conclusions"  are  suggested 

(1)  The  natural  directional  flux  density  distributions  J (0 * )  vs  0' 

(Pass  414,  Fig.  4.3  lower  right;  Pass  454,  Pig.  4.5  right)  are  roughly 
constant  in  space  and  time  in  good  agreement  with  early  published  results 
(Veber,  et  al  1962).  The  essential  characteristics  for  the  naturally 
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trapped  radiation  are 


Pass 

J(0' 

)  maximum 

.  Half-width  at 
ha If -maximum 
(h  w  h  m) 

4l4 

7 

x  107 

-2  -1  ,  -1 
cm  sec  ster 

21° 

454 

5 

x  10? 

-2  -1  .  -1 
cm  sec  ster 

28° 

The  maximum  flux  densities  J(0‘ )  maximum  noted  are  in  excess  of  those 
previously  reported  for  nearby  regions  of  space  because  of  corrections 
resulting  from  the  more  recent  calibrations  (Baicy,  et  al,  1963). 

(2)  The  ARGUS  directional  flux  density  distributions  for  the  trapped 
particles  are  more  variable  in  shape  both  with  time  and  space  location 
than  are  the  corresponding  natural  radiations.  The  ARGUS  J(6')  vs  6' 
curves  (Pass  414,  Pig.  4.4  lower;  Pass  454,  Pig.  4.6  right)  are  gener¬ 
ally  flatter  (near  the  maximum  of  the  curves)  than  for  the  corresponding 
natural  Radiation.  Especially,  the  J  (flO  vs  0*  half-width  at  half-maxi¬ 
mum  appears  to  be  a  variable  which  may  change  for  each  modulation  of 
the  count-rate  (see  CQba  vs  time  curves.  Pigs.  2.9,  2.10,  2.11  Lower). 

As  discussed  in  Sec  3*4  the  Important  factor  in  the  J(0')  vs  0'  shape 
variation  appears  to  be  the  variation  in  the  omnidirectional  penetrating 
background  radiation  which  is  shown  by  curves  B  in  Pigs.  2.10  and  2.11 
Lower.  The  essential  characteristics  for  the  ARGUS  trapped  radiation 
are 


Pass 

J(0')  maximum 

Half-width  at 
half-maximum 
(h  w  h  m) 

414 

(1.0-1. 9)  x  10®  cm^sec^ster’1 

24-30° 

454 

0.65  x  10®  cm“2sec"1ster‘1 

16° 
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It  is  necessary  clearly  that  a  great  deal  more  data  analysis  for 
ARGUS  trapping  is  required  to  understand  more  fully  the  directional 
flux  density  distributions  J(0')  vs  O'.  From  such  results  the  ARGUS 
mirror  point  distributions  would  be  obtained  immediately. 
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Fig.  4. 


L.  Geometry  involving  scintillation  counter  axis 

A 

direction  (  p.  unit  vector),  magnetic  flux 

*  A 

density  direction  (  Bj  unit  vector),  particle 
flux  density  direction  (specified  by  solid  angle 
3LQ  orientation)  and  various  angles. 


Co5a(  *  Coss'Stdo^  +*  CosQ 

and  aA'=  Cose'as'a^' . 
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«sr  . 


Pig.  4.2.  Response  of  S(a)  of  Detector  A,  directional  scintillation 


counter,  vs  the  angle  between  the  counter  axis  and  the 
direction  of  the  incident  p-radiation  (Baicy  et  al,  1962). 
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Pig.  4.3  Directional  count-rates  and  flux  densities  of  naturally 
trapped  radiation  vs  angle  6  (between  counter  axis  and  plane 
perpendicular  to  geomagnetic  field  direction)  obtained  by 
analysis  of  the  data  from  the  scintillation  counter  (Detector  A) 
on  Explorer  IV  a  few  minutes  before  entering  the  ARGUS  region 
(Event  I),  Pass  4l4. 

The  space  location  data  are  (for  250. 0  UT):  geographic  lati¬ 
tude  and  longitude,  3.8°  N  and  276°  E  respectively;  L,  1.40; 

B,  0 . U64  gauss;  altitude,  1781  km. 

Upper  and  lower  left.  Observed  count-rates  Co^g(8)  for  time 
intervals  noted  on  plots.  Count-rates  obtained  as  explained  in 
Chap.  2.  The  plotted  points  are  located  by  the  numbers  along 
the  curve;  numerals  1,  2,  3,  4,  5  indicate  count-rate  vs  6  for 
individual  modulations  as  shown  by  Pig.  2.9.  The  dashed  portions 
of  the  curves  are  theoretical  extrapolations. 

Upper  and  lower  right.  True  count -rate  C(6)  obtained  by  cor¬ 
recting  c0b8(®)  for  deadtime,  effective  area  and  counting  effi¬ 
ciency  of  the  scintillation  detector.  Directional  flux  densities 
J(0>)  obtained  as  explained  in  Chap.  4. 

The  J(0')  cutoff  and  half -width  at  half -maximum  (h  w  h  m)  are 
about  72°  and  21°  respectively  for  the  lower  right  plot. 
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Pig.  4.4  Directional  count-rates  and  flux  densities  vs  6  for  the 
scintillation  counter  during  penetration  of  the  ARGUS  region 
(Event  I)  a  few  minutes  later  than  for  Pig.  4.3.  C  (0), 

0  DS 

C(e),  and  J(0')  were  obtained  as  for  Pig.  4.3. 

The  space  location  data  are  (for  256.3  UT):  geographic  lati¬ 
tude  and  longitude,  17.6°  N  and  285. 5°  E  respectively;  L,  1.68; 

B,  0.221  gauss;  altitude,  1518  km. 

The  lower  case  letters  b,  c,  e  relate  the  data  of  Pig.  4.4 
with  that  of  Pig.  2.10  Lower  where  the  same  letters  are  em¬ 
ployed  for  corresponding  data. 

The  dotted  portions  of  the  curves  are  theoretical  extrapola¬ 
tions. 

The  J(0' )  cutoff  and  h  w  h  m  are  about  75°  and  24°  respectively 
for  the  lower  right  plot. 
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Pig.  4.5  Directional  count-rates  and  flux  densities  of  naturally 
trapped  radiation  vs  6  obtained  by  analysis  of  data  from  the 
scintillation  counter  (Detector  A)  a  few  minutes  before 
Explorer  I V  entered  the  ARGUS  region  (Event  II),  Pass  454. 
Cobs(e),  C(6),  and  J (©' )  were  obtained  as  for  Pig.  4.3.  The 
dotted  portion  of  the  curve  is  a  theoretical  extrapolation. 

The  space  location  data  are  (for  299* UT ) :  geographic  lati¬ 
tude  and  longitude,  4°  N  and  250°  E  respectively;  L,  1.35; 

B,  0.i57  gauss;  altitude,  1875  km. 

The  J (e* )  cutoff  and  h  w  h  m  are  about  60°  and  28°  respec¬ 
tively. 
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Pig.  4.6  Directional  count-rates  and  flux  densities  vs  0  for  the 
scintillation  counter  during  penetration  of  the  ARGUS  region 
(Event  II)  a  few  minutes  later  than  for  Pig.  4.5.  cobs(®)* 
C(0),  and  J(0')  were  obtained  as  for  Pig.  4.3.  The  dotted 
portion  of  the  curve  is  a  theoretical  extrapolation. 

The  space  location  data  are  (for  310,  UT):  geographic  lati¬ 
tude  and  longitude,  27.8°N  and  2J0°  E;  L,  2.08;  B,  0.269  gauss; 
altitude,  1405  km. 

The  J(0')  cutoff  and  h  w  h  m  are  about  40°  and  l6°  respec¬ 
tively. 
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Chapter  5 


Omnidirectional  Count-Rate  versus  Time.  Decays  of  Omnidirectional 
Radiation. 


5.1  Analysis  of  omnidirectional  count-rates,  Detectors  D  and  C, 
Explorer  IV 

An  analysis  was  made  of  the  data  obtained  through  Channels  1 

2 

(Detector  D)  and  3  (Detector  C),  the  shielded  (1.6  g/cm  ,  lead)  and 
unshielded  Anton  302  G-M  counters  respectively  (see  Pig.  2.1  block 
diagram).  Detector  C  had  no  intentional  shielding  beyond  that  of 
the  satellite  shell  and  contents  and  was  sensitive  to  electrons  of 
energy  greater  than  3  MeV.  The  shielded  Detector  D  had  an  electron 
energy  threshold  of  6  MeV. 

An  automatic  digitizer  was  used  to  analyse  the  telemetry  playback 
data  and  transform  it  to  punch  cards  and  a  computer  program  was  used 
to  obtain  true  count-rate  by  correcting  for  the  counter  deadtime, 

62.5  microseconds.  The  Explorer  IV  Passes  selected  for  analysis  were 
those  over  the  Huntsville  ground  station,  as  follows 

ARGUS  Event  I 


Pass 

Pass 

Pass 

August  27,  1938 

414 

415 

416 

August  28,  1998 

427 

August  29.  1998 

440 

August  30,  1938 

453 

454 

August  31,  1938 

466 
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September  1,  1958 
September  2,  1958 
September  3,  1958 
September  4,  1958 

ARCUS  Event  II 

August  30,  1958 
August  31,  1958 
September  1,  1958 
September  2,  1958 


479 

492 

505 

518 

453  454  455 

466 

479 

492 


(The  Passes  tabulated  above  are  arranged  In  3  columns  each  of 

which  corresponds  to  about  the  same  region  in  space  relative 

to  the  earth  -  the  Explorer  IV  satellite  orbital  frequency 
was  about  13  orblts/day). 

Pigures  5.1  and  5*2  are  sets  of  omnidirectional  count-rate  data 
plots  corresponding  respectively  to  ARQUS  Events  I  and  II.  The 
figure  captions  contain  explanatory  information.  The  curves  shown 
are  smoothed  averages  of  true  (corrected)  count-rates  with  the  actual 
points  not  Included. 

Pigures  5*3,  5*4  and  5*5  show  the  decay  of  the  ARQUS  shell  versus 
time  elapsed  T  since  Event  I,  or  II,  for  various  L  (B-L  geomagnetic 

coordinates)  values.  In  Pig.  5.3  the  data  points  for  Passes  505  and 

518  exhibit  increases  In  ARQUS  shell  intensity  (count-rates)  which 
are  due  probably  to  a  geomagnetic  disturbance  (Report,  State  University 
of  Iowa,  1959)* 


D,  Difference  |A  (unshielded)  -  B  (shielded)J  count-rates 
vs  time.  Assuming  that  electrons  only  are  counted  the 
difference  A-B  measures  roughly  the  omnidirectional 
counts/sec  due  to  all  electrons  of  energy  between  3  and 
6  MeV.  The  differences  A-B  are  for  count-rates  above 
background. 

T  ,  time  elapsed  (approximately)  since  Event  I. 
e 
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A-B  DIFFERENCE  10_2(counts/sec)  A/B  RATIO  10  (counts/sec) 
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Pigs.  5-2  True  count-rate  vs  time  for  omnidirectional  G-M  counters 
during  ARGUS  shell  penetrations  following  Event  II  for  Passes 
453,  454,  455,  466  and  479  of  Explorer  IV  Satellite  over 
Huntsville,  Alabama.  The  data  have  been  smoothed  so  that  the 
various  curves  are  averages  of  points  which  are  not  plotted. 

The  set  of  4  plots  for  each  Pass  are  arranged  as  follows: 

A,  unshielded  G-M  counter  (Detector  C,  Channel  3) 
which  counts  electrons  of  energy  greater  than  3  MeV; 
the  background  has  been  drawn  in  by  estimation. 

B,  shielded  G-M  counter  (Detector  D,  Channel  1) 
which  counts  electrons  of  energy  greater  than  6  MeV; 
the  background  has  been  drawn  in  by  estimation. 

C,  Ratio  of  A  (unshielded )/B  (shielded)  count-rates  vs 
time.  Assuming  that  electrons  only  are  counted  the  A/B 
ratio  can  be  expressed 


Ratio  A/B 


•i 

N(E)  dE 

3MeV 


where  N(E)  is  the  true  omnidirectional  counts/ sec 
per  energy  interval  in  MeV,  and  t  is  the  time  as  noted 
by  the  abscissa.  The  ratios  are  of  count-rates  above 
background. 
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D,  Difference  £a  (unshielded)  -  B  (shielded)]!  count-rates 
vs  time.  Assuming  that  electrons  only  are  counted  the 
difference  A-B  measures  roughly  the  omnidirectional 
counts/sec  due  to  all  electrons  of  energy  between  3  and 
6  MeV.  The  differences  A-B  are  for  count-rates  above 
background. 

T  ,  time  elapsed  (approximately)  since  Event  II. 

c 
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Pig-  5-3  Decay  of  ARGUS  shell  vs  time  elapsed  Tg  since  Event  I 
for  various  L  (Mclllwain  B-L  geomagnetic  coordinates)  values. 

The  plotted  points  are  true  (corrected)  omnidirectional 
count-rates  of  the  unshielded  G-M  counter  (Detector  C,  Channel  3) 
during  ARGUS  shell  Event  I  penetrations  of  several  Explorer  IV 
Passes.  The  count-rates  are  differences  between  total  and  back¬ 
ground  count-rates  due  to  electrons  (assuming  p-particle  trap¬ 
ping)  of  energy  greater  than  3  MeV. 

Data  for  Passes  492  and  505  were  taken  from  (Report,  State 
University  of  Iowa,  1959)  L  values  were  obtained  of  GE  TEMPO, 
Santa  Barbara,  California  through  the  offices  of  Major  Charles 
V.  Hulbert,  Defense  Atomic  Support  Agency. 
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Fig.  5.4  Decay  of  ARGUS  shell  vs  time  elapsed  Tg  since  Event  I 
for  the  Ratio  of  A  (unshielded  G-M  counter)/B  (shielded  G-M 
counter)  omnidirectional  count-rates,  for  various  L  values 
and  Passes  of  Explorer  TV  during  ARGUS  shell  Event  I  penetra¬ 
tions. 

The  Ratios  are  of  total  count-rates  minus  background.  For 
Pass  479  (not  shown,  but  refer  Fig.  5*3)  no  counts  above  back¬ 
ground  were  observed. 
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Pig.  5*5  Decay  of  ARGUS  shell  vs  time  elapsed  Tg  since  Event  II 
for  various  L  values.  The  plotted  points  are  true  (corrected) 
omnidirectional  count-rates  of  the  unshielded  G-M  counter 
(Detector  C,  Channel  3)  during  ARGUS  shell  Event  II  penetra¬ 
tions  of  several  Explorer  IV  Passes.  The  count -rates  are 
differences  between  total  and  background  count-rates  due  to 
electrons  (assumed)  of  energy  greater  than  3  MeV. 
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APPENDIX  1 


Power  Spectrum 


The  present  appendix  includes  an  analysis  of  the  roll  or  spin 
period  "anomaly"  discussed  briefly  in  Sec.  3.3  and  noted  by  Fig.  3.3. 

J.  Bock  (private  communication,  1964)  suggested  that  the  anomalous 
decrease  followed  by  increase  in  the  roll  period  observed  by  Naumann 
(1961)  might  be  the  result  of  mistaking  beat  rates,  between  roll  and 
tumble  motion  such  as  I  ^)*  2^  '  for  roll  rate  \jj  .  To 
evaluate  this  hypothesis,  a  trigonometric  series  analysis  of  the  Detector  A 
count-rate  vs  time  was  carried  out  for  the  time  interval  of  interest, 

Passes  4l4,  427  and  454.  Basically  the  theory  involved  is  as  follows 
(Blackman,  195") • 

If  a  function  X(t)  is  doflnc  1  for  all  times  t  in  the  interval  T  in 
which  periodicities  are  to  bo  Investigated,  the  power  spectrum  function, 

P  (  f)  can  be  written 


IS 


T/2 


X(‘)  e"1*™  dt 


P(f)»  l/T 

limT>  o 0  -j/2 


Ti ..  TT  ^  '•  ■’  '.’hi  ■  li  1  hi  of  P(f)  occur  yield  the  periodici* 

1  l>  .  I  x(i  ),  ^  o  moih.J  of  calculating  P (f )  involves  the 

mi  .. . .m*1  in  *  defined  as 

T/2 

C  (r)  *  li m 1/lT X (t ) X  (t  ^T)dT7 

T  -T/2 


(Al.l) 


(A1.2) 


where  T  is  called  the  lag  time.  Also,  P ( f )  is  the  Fourier  transform 
of  C(  T  )  and  so  may  be  written 


>(f) 


I 


-  i  2  TTi T 

C( t) e  dr  # 

QO 


For  functions  X(t)  which  are  defined  for  a  finite  time  interval 
but  not  for  all  t  (for  example,  the  Detector  A  count  rate  function),  an 
apparent  autocovariance  function  CQ(  T  )  may  be  defined  as 


C0  (T). 


and 


l/(T-  It 

fn 


r  (T-lrl)/e 

■ 

*  -(t-ItU/z 

m 


X  (t-r/2)  X (t  +  t/2) dt 
ItI  <  T  , 


C  (r)  *0;  for  Irl  >  T  I 
o  m  9 

where  <  Tn  (typically  Tm  is  Tn/'10).  Considering  Eq.  (A1.4), 
Eq.  (A1.3)  can  be  expressed  as 


P0(f) 


-12»f  T 

e 


dr  . 


(A1.3) 


(A1.4) 


(A1.5) 
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Equation  (A1.3)  applied  .to  a  very  long  record  containing  sinusoidal 
variations  with  time  of  discrete  frequencies  f^  (k  =  1,2,...)  would 
yield  P(f)  as  a  series  of  delta-function  spikes  at  these  discrete 
frequencies.  The  effect  of  using  records  of  finite  length  is  to 
broaden  these  sharp  spikes,  so  that  they  appear  in  the  form  (sin  x)/x 
where 


x 


2T<f-VTm 


The  function  (sin  x)/x  is  symmetric  about  x  =  0  with  maxima  at  x  =  0, 
7.725,  14.066,  20.371.  •••  of  corresponding  amplitudes  1,  0.1284, 
0.0709,  0.0490,  _ 

In  the  present  case  the  integrals  were  approximated  by  finite  sums 
with 


n-r 

C  *  l/(n-r>2l 

t-0 


X‘  X'  ;  r* 0,1,. ..m 

f  f+r 


where 

V  r  At, 

TmmAT> 

r-nAr, 

and 


x'*x-x,  "here  X  *X,(*At), 

t  f  1  1 


(Al. 6 ) 


(A1.7) 


(A1.8) 
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,1 

1' 


in  which  X  is  the  average  count-rate  defined  as 

n 

X-  l/n  2  X  .  1*1.9) 

H 

Subtraction  of  X  in  Eq.  (A1.8)  has  the  effect  of  reducing  the  zero 
frequency  peak  in  the  power  spectrum. 

To  interpret  the  power  spectrum,  it  is  noted  that  Eq.  (3.11)  namely 


6. 


ore  tin 


[belong  t+8)-  b*» M^t+/)coe(^  t+8)+&  o»s*£»+y)cos(^t48)] 

2  3 


(3.11) 


combined  wi£h  the  empirical  count-rate  function  used  by  Bock  (op.  clt. 
1964) 


(C  )  cos40 

Ob  9  0 


(A1.10) 


yield  an  empirical  analytic  expression  for  CQb8  vs  time  (see  Figs.  2.9 
and  2.10). 

Thus  CQto8  vs  t  (Pigs.  2.9  and  2.10  Lower)  can  be  expressed  in  a 
linear  series  of  trigonometric  functions  and  the  squares  of  the  ampli¬ 
tudes  of  these  terms  should  be  directly  proportional  to  the  Intensities 
for  the  frequencies  observed  in  the  directional  count-rate  data  and  also 
to  the  intensities  resulting  from  the  power  spectrum  analysis  assuming 

a  cos* 8  distribution  of  particle  flux.  Twenty  frequencies,  each  a 

•  • 

linear  combination  of  the  tumble  rate  and  roll  or  spin  rate  ^ 
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I' 


and  whose  amplitudes  are  not  zero  were  obtained  from  the  expan- 

4  -  *  •  •  • 

sion  of  cos  &.  The  rates  \j/,  <#>,3^and  3<£  have  zero 

amplitudes  and  thus  should  not  appear  in  the  power  spectrum.  The  in¬ 
tensities  corresponding  to  the  remaining  frequencies  are  a  function  of 

A  A 

b^  =B*Lonly  (see  Sec.  3.1).  These  frequencies  and  relative  inten¬ 
sities  are  listed  in  Table  Al. 


Table  Al  shows  that  dominant  frequencies  are  2  4>,  Zt, 

i  •  #  « 

4>  +  2  s  f  l<^>-2^l  s  f  •  The  relative  intensities  corres¬ 
ponding  to  these  frequencies  are  plotted  as  a  function  of  b^  in  Pig.  Al. 
Figure  Al  shows  that  when  b.  =  0  forCOS(B*L)  s  9  0°1  both 
the  2  (p  and  2  ^  terms  have  equal  intensity  while  the  f  and  f _ 

A  A 

terms  have  zero  intensity.  As  the  angle  between  L  and  B  decreases  from 

90°  to  0°  or  increases  from  90°  to  180°  the  Intensity  for  2  (f>  decreases 

A  A 

monotonically  to  zero  which  is  to  be  expected  when  OTCOS(B*L) 
approaches  0°  or  180°  in  which  cases  b^  and  b^  go  to  zero  and  b^  goes  to 
unity.  By  Eq.  (3.11),  if  bg  *  b^  =  0,  Q  is  Independent  of  <f>  and 

b^  =  |  and  physically  the  count-rate  would  show  no  variation  due  to 

A  2 

rotation  about  L.  It  is  also  clear  by  Table  Al  that  for  b^  =  1  all 

#  •  0 

intensities  involving  (p  go  to  zero  and  2  ^  and  alone 

are  present. 

* 

The  behavior  of  the  intensity  for  frequency  2<£  was  thus  expected 
physically.  However  the  variation  of  the  intensities  of  2  \Jf  % 
and  f  with  b^  could  not  be  anticipated  from  physical  considerations. 

It  is  seen  (Pig.  Al)  that  the  intensity  for  2</r  decreases  to  zero 
for  a  value  of  about  . 553  for  b^ ;  this  corresponds  to  arcos(B#L) 
about  56.5°.  The  intensity  for  2\Jf  then  rises  to  a  very  large  value 
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is  almost  the  reverse.  They  are  zero 


about  49°.  Since  this  latter  intensity  is 


for  b^  =  1.  The  behavior  of  the  intensities  for  f  .  =  <£*2  * 

and  f _  =  If- 2  *  I  is  almost  the  reverse.  They  are  zero 

at  b^  =  0  and  1  and  are  a  maximum  for  b^  about  .668  corresponding  to 
OrCOS(B  •  L)  about  49°.  Since  this  latter  intensity  is 

near  its  maximum  when  the  intensity  for  2  ^  is  zero,  it  is  seen 

* 

that  the  frequencies  f+  and  f_  are  dominant  over  2\fr  in  the  range 
.33  <  b^  «  .79;  while  for  b^  outside  this  range,  2^  is 
dominating  frequency. 

For  Pass  454  during  the  time  interval  from  300.0-305.0  UT  b^  was 
calculated  to  vary  from  .01  to  0.36.  It  was  during  this  interval  that 

the  long  term  periodicities  mentioned  by  Naumann  (1961)  representing 

•  • 

2*  were  measured.  Figure  A1  shows  that  2+  should  be  the  domi¬ 
nating  frequency.  However  for  Pass  4l4,  b^^  was  calculated  to  vary 
from  0.33  to  O.70  in  the  interval  from  250.0-255*0  UT,  Fig.  A1  shows 
that  1^  ±  2  \j/ 1  would  be  the  dominating  frequency.  In 

Pig*  3*3  the  frequency  l<£-2  \£l  is  plotted  with  squares  for 
Pass  414,  427,  and  453. 

The  appearance  of  harmonic  or  beat  frequencies  in  the  count-rate 
vs  time  data  appears  to  account  for  the  apparent  anomalous  behavior 
of  the  roll  period  vs  time  as  determined  by  Naumann  and  plotted  as 
circle  points  in  Fig.  3*3  (in  the  interval  August  24-29  approximately). 
J.  Bock  (private  communication,  1964)  identified  the  dominant  frequency 
in  the  anomalous  region  as  l$-  2*1  rather  than  2* 

and  the  present  analysis  agrees  with  his  Identification. 

An  example  of  the  power  spectrum  analysis  is  given  by  Fig.  A2  for 
Pass  4l4  which  falls  within  the  region  of  apparent  anomalous  behavior. 
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The  frequencies  and  relative  intensities  of  the  dominant  peaks  are 
peak  frequency  (cycles  sec'1)  relative  intensity 


1 

0.0080 

703 

2 

0.0430 

618 

3 

0.0628 

141 

4  ‘ 

0.1859 

108 

5 

0.2935 

138 

6 

0.3306 

296 

7 

0.4792 

114 

Peak  1  is  not,,  physical  but  results  from  the  finite  range  of  integration 

over  T  .  Side' lobes  of  peak' 1  (due  to  the  sin  x  behavior  discussed 

x 

above)  appear  to  have  interfered  with  peak  2  increasing  its  intensity 

somewhat.  Peak  3  is  interpreted  as  a  side  lobe  of  peak  2  since  a  side 

lobe  calculated  by  Eq.  (Ai»6)  occurs  for  f  =  0.0635.  Peak  5  is  at  the 

» 

frequency  corresponding  to  2*  the  tumble  period  being  about  7 

9 

seconds.  Since  peak  6  appears  at  a  larger  frequency  than  for  2  4>  and 

has  a  greater  intensity  than  for  2*  it  must  be  interpreted  as 

•  »  * 

associated  with  the  frequency  <^>+2  because  2^  never 

exceeds  2  for  Explorer  IV  after  the  propeller-like  motion  is  es¬ 


tablished,  and  from  Pig.  A1  no  other  frequency  could  exceed  that  for 
2  <t>  in  intensity  except  at  large  values  of  b1  when  the  2^>  Peak 
would  be  buried  in  noise.  Peak  6  could  not  be  associated  with 

1  ^  I  *  •  • 

since  if  yjf  « <f)  one  can  show  |^>—  2^l«  <j>  The 

interpretations  of  peaks  5  and  6  of  Pig.  A2  yield 

tumble  frequency  *  0.1468  cycles  sec"1 

\ff  roll  (or  spin)  frequency  *  O.0919  cycles  sec"1  (Al.ll) 

Prom  Eq.  (Al-11),  as  a  cross  check 


frequency  corresponding  to  =  0.0371  cycles/sec 

frequency  corresponding  to  2\j/  =  0.1838  cycles  /sec 

frequency  corresponding  to  2^+2\fr  *  0.4774  cycles/ tec  (Ai.12) 
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The  assignments  above  allow  the  identification  of  peak  2  (.0430  cycles 

sec"'*')  as  due  to  I* -2*1,  peak  4  (.1859  cycles  sec"*)  as 

«  '  *  • 

due  to  2  +  and  peak  7  (-4792  cycles  sec  * )  as  due  to  2**2* 

The  values  of  Eq.  (Al.ll)  compare  with  Bock's  more  accurate  values  for 
Pass  4i4  as  follows 

(j)  tumble  frequency  =  0.14490  cycles  sec”* 

roll  (or  spin)  frequency  =  0.09382  cycles  sec”*  (A1.13) 

Prom  Eqs.  (Al.ll)  and  (A1.13),  it  is  seen  that  about  a  2  percent  differ¬ 
ence  exists..  Using  Bock's  values,  Eq.(A1.13),  yields  for  frequency 

«  j* 

corresponding  to  2^*  r  =  0.0427,  in  much  better  agreement 

with  the  peak  2  value  than  as  stated  by  Eq.  (A1.12).  However  the  fre¬ 
quencies  obtained  by  a  power  spectrum  analysis  would  be  close  enough  to 
the  correct  values  to  be  useful  as  starting  parameters  in  the  "closure 
search"  procedure  described  in  Sec.  3.3. 

It  appears  possible  also  to  perform  an  orientation  analysis  with  the 
power  spectrum.  If  a  separate  power  spectrum  were  computed  for  each 
minute  of  a  pass,  and  the  relative  intensities  of  the  relevant  peaks 
examined  by  use  of  Pig.  Al,  an  average  b^  for  that  minute  could  be  ob¬ 
tained.  A  plot  of  b^  vs  time  then  could  be  made  and  a  determination  of 

A 

the  absolute  direction  of  L  carried  out.  Such  procedure  would  be  the 
same  type  as  done  by  Bock  (1964)  except  that  he  used  an  analog  device 
in  which  only  count  rate  maxima  were  involved.  The  method  here  described 
would  use  all  the  data  and  might  possibly  yield  better  results,  that  is 
if  the  signal-to-noise  ratio  were  large  enough.  Figure  A2  is  the  power 
spectrum  for  the  entire  Pass  4l4  consequently  the  intensities  have  been 
averaged  over  about  9  minutes,  during  which  b^  varies  from  0.3  to  about 
0.8. 
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Table  Al.  Harmonics  and  beat  frequencies  of  tumble  rate  <^>  and  roll  (or  spin) 

*  4 

rate  ^  obtained  by  expansion  of  cos  0  (Eq.  (3-11))  in  a  Fourier  type  series. 


Harmonics  and 
beat  frequencies 

Calculated  intensity 

2  4 

V 

.  1/28 

(1-b2)2  (5-3b2)2 

4  4 

V 

-  9/212 

(1-b2)4 

CM 

*3 

•  1/28 

5-l8b^  +  5bJ  2 

4  *f/ 

•  V 

-  !/212 

2  4 

3-30b*  +  35b^ 

4>+2\j/  -f+ 

V 

.1/2* 

b2  (1-b2)  (5-b2); 

1  4  -  2  +  1  *  f_ 

*6  ; 

=  I5 

•  l 

<(>+4*1/ 

*7  ' 

.1/26 

b2  (l-b*)(3-7b*): 

14-4+1 

I8 

-1? 

2  4  *2  + 

V 

■  1/28 

(1-b2 ) 2  (1+b2)2 

\z4~  2.+\ 

I10 

•!9 

2+*4+ 

*11 

-  1/210  (1-b?)2  (1-7D? 

I  z4  -  4+\ 

I 

a  I 

•  f 

A12 

A11 

3<*>‘  +  2</» 

J13 

*  1/26 

b2  d-b2)5 

l3<#>-  2*1 
•  • 

*14 

■  X13 

3+ +4^ 

I15 

-1/28 

u-b2)3 

13  *-4  + 1 

*16 

•hi 

44*2+ 

I17 

-  1/210  <l-b2|4 

\4+-2+\ 

I18 

~X17 

4<j>+4  \j/ 

J19 

-  1/214  (1-b2)4 

14+-4+1 

*20 

■X19 
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Pig.  Al.  Several  harmonic  and  beat  frequency  intensities,  from 
Table’  Al,  vs  b^.  The  frequencies  involved  are  A,  2^,. 


B£j>;  DjiiZtii  C,l2<£i2^l. 

The  ordinates  have  been  multiplied  by  2 


12 
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INTENSITY  IN  ARBITRARY  UNITS 


A.  2  Power  Spectrum  Analysis  for  Pass  4l4.  Relative  power 
or  intensity  P(f)  is  plotted  vs  tumble  and  roll  (or  spin) 
beat  frequencies  or  harmonics.  The  numbered  peaks,  1  to  7» 
are  explained  in  the  text. 


FREQUENCY  (CYCLES  SEC 


Distribution  List 


1.  Director,  U.  S.  Army  Ballistic  Research  Laboratories,  ATTN:  Mr.  Ed  Baicy, 

Aberdeen  Proving  Ground,  Maryland  21005 

2.  Commanding  Officer,  U.  S.  Army  Research  Office  -  Durham,  Box  CM,  Duke 

Station,  ATTN:  Dr.  R.  Mace,  Durham,  North  Carolina  27706 

3.  AFWL  ATTN:  WLRP;  WLRB;  WLAX,  Kirtland  APB,  New  Mexico  87117 

4.  AFCRL  ATTN:  Dr.  L.  Katz,  L  G  Hanscom  Pld. ,  Bedford,  Mass.  01731 

5.  Director,  Goddard  Space  Plight  Center,  National  Aeronautics  and  Space 

Administration,  ATTN:  Dr.  V.  Hess  (3  cys),  Greenbelt,  Maryland 

5.  Director,  Advanced  Research  Projects  Agency,  Washington,  D.  C.  20301 

(5  cys) 

6.  DDC  Cameron  Stn. ,  Alexandria,  Virginia  (20  cys) 

7.  Director,  Defense  Atomic  Support  Agency,  ATTN:  Document  Library  (5  cys); 

RAEL,  Washington,  D.  C.  20301 

8.  RAND  Corporation,  1700  Main  Street,  Santa  Monica,  California,  ATTN: 

Dr.  P.  Tamarkln 

9.  Aerospace  Corporation,  P.  0.  Box  95085,  Dos  Angeles,  California,  ATTN: 

Dr.  S.  Preden,  Dr.  J.  Vette,  Dr.  R.  Harrison 

10.  Lockheed  Missiles  and  Space  Company,  A  Division  of  Lockheed  Aircraft 

Corporation,  1111  Lockheed  Way,  Sunnyvale,  California,  ATTN:  Dr.  M. 
Walt,  Dr.  R.  Smith,  Dr.  J.  Caldls,  Dr.  W.  Iahof,  Dr.  R.  Meyerott 

11.  Lawrence  Radiation  Laboratory,  P.  0.  Box  808,  Livermore,  California,  ATTN: 

Dr.  C.  Mclllwaln,  Dr.  H.  West,  N.  Christofllos,  Dr.  G.  Bing 

12.  State  University  of  Iowa,  Administration  Building,  Iowa  City,  Iowa,  ATTN: 

Document  Control  Station  for  Dr.  J.  Van  Allen 

13.  Bell  Telephone  Laboratories,  Whlppany  Road,  Whippany,  New  Jersey,  ATTN: 

Dr.  W.  Brown,  Dr.  C.  Roberts 

14.  North  American  Aviation,  Space  and  Information  Systems  Division,  12214 

Lakewood  Boulevard,  Downey,  California 

15.  The  Catholic  University  of  America,  Washington,  D.  C.  20017,  ATTN: 

Dr.  C.  Chang 


Distribution  List  (continued) 


16.  General  Electric  Company,  Technical  Military  Planning  Operation,  735 

State  Street,  Santa  Barbara,  California,  ATTN:  Dr.  R.  Hendrick, 

DASA  Data  Center  (2  cys) 

17.  IIT  Research  Institute,  10  West  35th  Street,  Chicago,  Illinois  60616, 

ATTN:  Dr.  B.  McCormac,  Dr.  L.  Reiffel 

18.  Los  Alamos  Scientific  Laboratory,  P.  0.  Box  1663,  Los  Alamos,  New  Mexico, 

ATTN:  Dr.  A.  Petchek,  Report  Librarian 

19.  General  Atomic,  A  Division  of  General  Dynamics  Corporation,  P.  0.  Box  1111, 

San  Diego,  California  92112,  ATTN:  Dr.  David  B.  Chang 

20.  Director,  Applied  Physics  Laboratory,  John  Hopkins  University,  8621 

Georgia  Avenue,  Silver  Spring,  Maryland,  ATTN:  Dr.  R.  Pise he 11 

21.  Air  Force  Technical  Applications  Center,  Washington  25,  0.  C.,  ATTN: 

Mr.  D.  Northrup 


